Optical, vibrational, thermal, electrical, damage and phase-matching properties of 

lithium thioindate 



O 

o 

> 

CJ 

o 

00 

(N 



C/3 

o 

> 

q 
o 



(N 
> 
(N 
oo 
O 

o 

cn 

o 
o 

• 1— I 



X 



S. Fossier, S. Salaiin, J. Mangin and O. Bidault 

Laboratoire de Physique de I'Universite de Bourgogne (UMR-CNRS 5027), 9 avenue A. Savary, F-21078 Dijon, France. 

I. Thenot and J.-J. ZondjB 

BNM-SYRTE, Observatoire de Paris (UMR-CNRS 8630), 
61 avenue de I'Observatoire, F-75014 Paris, France. 

W. Chen 

Laboratoire de Physicochimie de I'Atmosphere (UMR-CNRS 8101), 
Universite du Littoral Cote d'Opale, 145 Av. Maurice Schumann, F-59140 Dunkerque, France. 

F. Rotcrmunc0 and V. Petrov 
Max-Bom Institute for Nonlinear Optics and Ultrafast Spectroscopy, 2A Max-Bom St, D-12489 Berlin, Germany. 

P. Petrov and J. Henningsen 
Danish Institute of Fundamental Metrology, B 307 Matematiktorvet, DK-2800, Lyngby, Denmark. 

A. Yelisseyev, L. Isaenko and S. Lobanov 
Design and Technological Institute for Monocrystals SB-RAS, 
43 Russkaya St, RU-630058, Novosibirsk, Russia. 

O. Balachninaitell G. SlekysU and V. SimtkaitifU 

Altechna Co. Ltd, Ukmerges 41A-203, LT-2004, Vilnius, Lithuania. 
(Dated: February 2, 2008) 

Lithium thioindate (LiInS2) is a new nonlinear chalcogenide biaxial material transparent from 0.4 
to 12 /^m, that has been successfully grown in large sizes and good optical quality. We report on 
new physical properties that are relevant for laser and nonlinear optics applications. With respect 
to AgGaS(e)2 ternary chalcopyrite materials, LiInS2 displays a nearly-isotropic thermal expansion 
behavior, a 5-times larger thermal conductivity associated with high optical damage thresholds, 
and an extremely low intensity-dependent absorption allowing direct high-power downconversion 
from the near-lR to the deep mid-IR. Continuous- wave difference- frequency generation (5-11 ^.m) of 
Tiisapphire laser sources is reported for the first time. 

PACS numbers: 42.70.Mp; 42.65. Ky, 42.65.Lm, 65.40.Ba, 65.40.De, 78.20.Fm, 78.20.Hp, 78.20.Jq, 78.20.Nv, 
78.30.Fs, 61.80.Ba 



I. INTRODUCTION 

The intensive search for new nonhncar materials for 
the generation of coherent tunable radiation in the mid- 
IR (2 — 20 /im), a spectral range of importance for molec- 
ular spectroscopy, atmospheric sensing and various opto- 
electronic devices, remains a continuing challenge. In 
contrast with the oxides (with IR transmission range 
not exceeding ^ 5 /im) that attracted much attention 
lately, such deep mid-IR materials cannot be grown by 
well mastered hydrothermal, flux or Czochralski meth- 
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ods. Instead the more complex Bridgman-Stockbarger 
growth technique in sealed (high-pressure) ampoules, 
with volatile and chemically reactive starting compo- 
nents, is the only method used to produce large size sin- 
gle domain crystals, and this certainly hampered their 
development all the more that special post-growth treat- 
ments are needed to restore stoichiometry and improve 
their optical quality. 

To date, only few suitable nonlinear crystals combin- 
ing an extended transparency from the visible/near-IR to 
the deep mid-IR and large enough birefringence to allow 
phase-matching over the whole transparency range are 
available. The majority of these compounds belong to the 
ternary chalcogenide semi-conductors of the A^B^^^CY' 
family, where A=Cu, Ag; B=A1, Ga, In; C=S, Se, Te |l[. 
Except for AgInS2, which can exist in both the tetrago- 
nal and orthorhombic phases these crystallize in the 
chalcopyrite (CuFeS2) structure with J42d space group 
and tetragonal symmetry (class) 42m but only few of 
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them possess sufBcient birefringence. Most of the com- 
merciahy available nonlinear crystals (AgGaS2 or AGS, 
AgGaSea or AGSe, ZnGePa or ZGP, GaSe, CdSe) or 
more confidential ones (HgGa2S4 or HGS, AgsAsSa or 
proustite, TkAsSes or TAS) have their advantages and 
drawbacks Despite their higher nonlinearity com- 
pared with the oxides, one limitation stems from their 
low energy bandgap, which prevents the use of pump 
lasers below 1 /im due to severe linear and two-photon 
absorption and associated thermal effects. ZGP has a 
high nonlinearity (72 pm/V) and thermal conductivity 
but requires pump wavelength above 2 /im due to its 
residual absorption in the near-IR. The noble metal com- 
pounds (AGS, AGSe) have the lowest residual absorp- 
tion, but their poor thermal conductivity due to low- 
energy phonon spectrum limits their performance, espe- 
cially in cw applications Furthermore, their 
thermal expansion is anisotropic along the direction par- 
allel and perpendicular to the optic axis, which is a source 
of thermo-mechanical stresses. CdSe's transparency ex- 
tends up to 18/Ltm but its low birefringence and nonlin- 
earity limit its phase-matching capability. GaSe is a soft, 
cleaving compound that cannot be cut at directions dif- 
ferent from the optical axis, further its large birefringence 
results in severe walk-off limitations. 

Two new materials belonging to the A^B^^^C^^ chalco- 
genide family, where the metal cations are replaced with 
the lighter alkali metal (A=Li), can now be added to 
this limited list: LiInS2 or lithium thioindate (LIS) and 
LiInSe2 or lithium selenoindate (LISe) . They were shown 
to crystallize with the /3-NaFe02 structure (orthorhom- 
bic Pna2i = Cf^, symmetry, point group mm2), which 
is a distorted superstructure of the wurtzite type, as 
early as 1965 0- Hence, contrary to the other mid-IR 
compounds, these new materials are biaxial and isostruc- 
tural to, e.g., KTiOP04 (KTP). The linear and nonhn- 
ear properties of LIS were briefly studied by Boyd and 
co-workers in the early 70's, using small and poor qual- 
ity samples |^] , but very few works (mainly focused on its 
optical ban dgap characterization and photoluminescence 
spectra [IS iH El IH Q 113 ) have followed since until 
recently due to problems associated to its reliable growth 
in larger sizes allowing to design good optical quality sin- 
gle domain elements for practical applications . Thin 
films of LiInS2 on glass and silicon substrates were also 
characterized 17j. The growth problems, associated with 
the high chemical activity of Li, were solved recently ^3 
and large crystals (up to 20 mm in length) of high op- 
tical quality became available for practical nonlinear op- 
tical applications 19,, JH, JlJ. To our knowledge, there 
is only one other group which has succeeded previo usly 
in synthesizing this material with comparable size |22j| . 
but this work was also stopped recently. While only non- 
phase-matched second-harmonic generation (SHG) was 
used by Boyd and co-workers to determine the nonlinear 
properties of LIS 2], the first phase-matched SHG ex- 
periment using the tunable output of a pulsed LiNbOa- 
OPO near 2.5 fj,m was reported only recently The 



range of fundamental wavelengths was extended to 6 /im 
using the radiation from a free-electron laser and ZGP- 
frequency doubler |23| and rough estimation of the ef- 
fective nonlinear coefficients was extracted from those 
measurements. Direct parametric down-conversion in the 
range 4.8—9 fim, using a femtosecond (200 fs) Ti:sapphire 
amplifier as the pump source, was simultaneously re- 
ported with an upper-bound value of the two- pho ton ab- 
sorption (TPA) at 800 nm of /3=0.04 cm /GW I2ll owing 
to the high bandgap energy of LIS |T1. Ill I2I l23t . Such 
a value is about 100 times less than the value for AGS or 
1000 times less than for AGSe at 1.32 /im p3 |. 

These preliminary experiments allowed to assess the 
potential of LIS over the existing mid-IR materials, and 
justify the need for a complete characterization of its op- 
tical properties. In this paper, we report new results 
from an intensive investigation campaign involving a net- 
work of laboratories in Europe, on the characterization of 
the main properties relevant for laser applications: struc- 
ture and transmission (Section^, specific heat, ther- 
mal expansion and thermo-optic coefficients fSection lIII|l . 
accurate lattice phonon spectra (Section II V|) . electro- 
optic and piezo-electric constants (Section 0, linear op- 
tical and thermo-optic dispersions fSection IVI|) . phase- 
matching predictions (Section lVII|) . nonlinear coefficients 
fSection I VIIip . parametric down-conversion fScction |lXp 
and optical damage thresholds (Sectional. The results 
presented here show that LIS is superior to the existing 
chalcopyrites in terms of high-power deep mid-IR down- 
conversion applications pumped directly from the near 
IR. The preliminary work on LISe has just started and 
will not be reported here [2B |. 



II. COMPOSITION, STRUCTURE AND 
TRANSMISSION 

Single-domain crystals of LIS are grown by the 
Bridgman-Stockbarger technique in a vertical setup with 
counter-pressure on seeds along (001) and (010) direc- 
tions H M 113 III El m m . The directed crystalfiza- 
tion is performed from a melt of elemental Li, In and 
S. In the melt zone the temperature is maintained at 
^ HOOK and decreases to ~ 900 K in the growth zone. 
The charged ampoules are moved in an optimal thermal 
gradient of 10 — 15 K/cm at a rate of ~ 10 mm/day. This 
technique allows the growth of crystal ingots with diame- 
ters up to 20 mm and lengths up to 50mm . Although 
sometimes colorless to a great extent, the as-grown in- 
gots are usually milky because of small inclusions of var- 
ious phases and it is necessary to anneal them in Li2S 
or S2 vapor at temperature close to the melting point 
{Tmelt ~ 1000°C) 18J. The crystal coloration after the 
thermal post-growth treatment changes then from almost 
colorless (or slightly yellow) for as-grown, to salmon-rose 
tinge depending on the growth and post-growth treat- 
ment conditions and chemical composition. Optical ele- 
ments with aperture up to 1 cm^ and length up to 15 mm 
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FIG. 1: Fragment of the orthorhombic unit cell structure 
(the box frame gives the orientation of the unit cell). The 3 
empty areas materialized with light grey spheres indicate the 
location of octapores. 



could be prepared for optical measurements or for non- 
linear frequency conversion. 

Precise chemical analysis - using differential dissolu- 
tion techniques (DDT) and inductively-coupled plasma 
(ICP) method - showed that the composition of color- 
less (or yellow tinge) and rose samples is Lip.cislni .n.'^ S2.n n 
(Lio.gelni.oeSi.gs) and Lio.86ln1.04S1.98 respectively Ts^. 
Lithium deficit (A^Li • -^In < 1) is observed for all sam- 
ples, the annealed rose ones being characterized by a 
higher sulfur content (A'^Li+in : -^S < !)■ The band- 
structure nature of LIS has been confirmed to be of 
the direct typ e (direct transitions between parabolic 
bands) [11 Il4j . The forbidden band gap values esti- 
mated by us [Eg = 3.72 and 3.57 eV at 80K and 300K, 
respectively, which are in good agreement with previous 
measurements ^1 0, 1(22, 23] ) , the position of the long- 
wave edge of transparency (13/im at an absorption of 6 
cm~^) and the crystal symmetry group were found iden- 
tical for all crystal colorations although small variations 
in latti ce p arameters due to point defects were also ob- 
served [23 ■ 

A detailed and general structural analysis of A^B^^^C^^ 
compounds, where A^ is an alkali metal such as Li, has 
been given by Kish et al [sj- A fragment of the struc- 
ture of LIS is shown in Fig. [l] (see Figs. 2.1.1-2.1.2 
in Ref. 28] for the full unit cell representation). It is 
formed by LiS4 and InS4 tetrahedrons and the S^^ ions 
are arranged in hexagonal packing with tetragonal and 
octahedral cavities (tetra- and octapores) [3ll. Com- 
pared to the denser chalcopyrite structure of AGS or 
AGSe [s^, the /3-NaFe02 structure of Li-chalcogenides 
allows the presence of empty cavities within the unit cell 
volume. The less dense /3-NaFe02 structure allows in 
particular doping with rare-earth (RE) active ions which 
can be hosted inside the octapores. Such a preliminary 
doping experiment with Nd'^"'" ions was previously re- 
ported with the aim of investigatingREiLIS potential 
as a self down-converting material 33]. The values of 



FIG. 2: Polarized transmission spectra of an annealed LIS 
sample (thickness 5 mm), not corrected for the Fresnel loss. 



the orthorhombic lattice parameters are a = 6.874(1)A, 
6 = 8.033(2)A, c = 6.462(1)A for as-grown colorless LIS, 
and a = 6.890(1)A, b = 8.053(1)A, c = 6.478(1)A for 
slightly yellowish samples. They are slightly lower than 
the values reported by Hoppe [1] and Boyd et al 0] or 
Kamijoh et al 0, and Kish et al [s^], but most 
closely to the last one. They modify to a = 6.896(1)A, 
6 = 8.058(2)A, c = 6.484(4)A for rose annealed samples. 
Concomitantly, the density changes from p — 3.52g/cm^ 
to p = 3.44g/cm'^. We have chosen c < a < h where 
c is the polar two-fold axis which coincides with the se- 
lection of Boyd et al This convention satisfies the 
recommendation of Roberts |35j for unique designation 
of the orthorhombic class mm2. Note that in the fol- 
lowing, in contrast to Roberts [s^ and the ANSI/IEEE 
Standard |36l ]. we do not introduce an additional frame 
for reporting the nonlinear tensor properties but for sim- 
plification and in accordance with the tradition use the 
crystallographic ahc frame for reporting the nonlinear op- 
tical susceptibility (see Section rVII|) . Also, instead of 
small case letters xyz, as suggested by Roberts [s^, we 
use the capital letters XYZ for designation of the prin- 
cipal optical axes such that nx < ny < nz- Since for 
LIS one has rif, < Ua < ric the principal frame assignment 
will be X <-> fe, y ^ a and Z c. 

Fig. 121 shows the polarized transmission spectra of a 
thick {L — 5 mm) annealed optical element. Depending 
on the growth and annealing conditions, some absorp- 
tion bands due to the various point defects - in partic- 
ular excess sulfur ions occupying interstitial sites - may 
limit the transmission window in the visible side. These 
visible bands, responsible for the rose coloration of the 
annealed samples, restrict the range of high-transmission 
(T ~ 70%) to 0.8 - S.lfim. The post-growth thermal 
treatment, however, improves considerably the trans- 
parency within the above window. The absorption co- 
efficient a{X) within the full transparency window de- 
duced from the spectra amounts to 0.1 - 0.15 cm~^. 
These values are overestimated, since from direct laser 
line transmission experiments using various laser sources 
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in the range 0.78 — 2.53 fj,m and several samples the aver- 
age residual absorption of a typical annealed sample was 
found to be lower. At 2.53 /im, e.g., the total loss coef- 
ficient (absorption and scattering) estimated amounts to 
0.05 cm^^. At 780 nm the loss coefficient measured with 
a diode laser source ranged between 0.05 and 0.2 cm^^ 
depending on the strength of the weak point-defects re- 
lated bands 0| seen also in Fig. |21 



III. THERMAL AND THERMO-OPTIC 
PROPERTIES 

The thermally-related properties of a new nonlinear 
material are crucial in assessing its potential in real non- 
linear conversion devices pumped by high power cw or 
pulsed lasers. The performance of devices based on mid- 
IR chalcogenides is often limited by deleterious thermal 
effects (lensing, expansion) ^ ■ The thermal fig- 

ure of merit of a material can be defined by the quantity 
J] = (dn/dT)/K where dn/dT is the thermo-optic in- 
dex variation and K is the thermal conductivity. The 
larger this quantity the more severe are the thermal lim- 
itations. A second thermal quantity, td = wgCp/AK 
where Cp is the mass specific heat of the material and wa 
is the beam radius, gives the characteristic time of heat 
diffusion outside the pumped volume. The larger this dif- 
fusion time, the stronger the thermal lensing effects 
It is hence important to know Cp, K and dn/dT. In addi- 
tion, the thermo-optic coefficients are needed to predict 
phase-matching performance in x^^^ interactions. Let us 
point out that no data has been published on the thermal 
properties of single domain crystalline LIS. 



A. Specific heat at constant pressure 

Using AI2O3 as a reference we estimated Cp{T — 
SOOAT) = 92.9 ± l.lJ/mol/K for crystalline rose colour 
LIS (57.91 mg) which is only slightly higher than what 
has been previously measured for polycrystalline LIS 
(90 J/mol/K) in Ref. [s^- The temperature dependence 
in the [-160°C; 210°C] range could be fitted by 



Cp(T) = 114.02 - 5708 x T^^ - 187540 x 



(1) 



where T is in Kelvin. The general theoretical expression 
of Cp (T) for anharmonic solids is 



Cp (T) = 12R 



N 



k=l 



(2) 



where R is the molar gas constant, xd = To/T {To is 
the Debye temperature) and F{xd) is the Debye function 
describing the temperature dependence in the harmonic 
lattice vibration approximation. 
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The value of N and the absolute magnitudes of the coef- 
ficients Cfc in the sum expansion in ^ can be considered 
as a measure of the degree of lattice anharmonicity, while 
the sign of Ck is essentially determined by the shape of the 
interatomic potential. In Eq. (0) , the leading coefficients 
ci and C3 were found negative for all A^B^^^C^^ chalco- 
genide compounds, and the magnitude of Ck was found 
^ 10 times higher and independent of the anion C^^ for 
Li-compounds. This means that the anharmonic contri- 
bution of the lattice potential energy is much stronger in 
LilnCyi compounds (iV = 4 in Eq.EJ) than in A^^'^^^G^^ 
chalcopyritc compounds {N = 3). In LIS, this is re- 
lated to the specific nature of the Li-S bond, which is 
about twice weaker than the In-S bond [s^- As a conse- 
quence, the temperature dependence of Cp(LIS) is much 
weaker than for AgInS2 and AGSe [39|, and the result- 
ing Cp{T) values for LIS are the lowest of all chalco- 
genides. By a calorimetric absorptance {A) measure- 
ment the absorption coefficient a can be estimated from 
A = mCpAT/PAt = 1 — exp(— aL) where m is the mass 
in mol, AT is the temperature rise, P is the laser power 
and At is the irradiation time. We irradiated at 1064 nm 
an annealed LIS sample of dimension 5x5x5 mm"^ cut at 
28° from the c crystallographic axis for light propagation 
in the a-c plane (mass m = 0.438 g = 2.356 x 10"'^ mol) 
and obtained an absorption coefficient of a = 0.037 cm~^ 
(compare Section llT)l . 



B. Thermal expansion 

The principal thermal expansion and thermo-optic co- 
efficients were determined by using the absolute interfer- 
ometric dilatometer and the experimental procedure de- 
scribed in Ref. [i^] . A two-beam modified Mach-Zehnder 
interferometer arrangement is employed for dilatation 
measurements, while thermo-optic coefficients are ob- 
tained using the same set-up by recording changes in op- 
tical thickness of a sample acting as a thermal scanning 
Fabry-Perot interferometer, with the natural refiectivity 
of the two opposite parallel facets. Three parallelepipedic 
samples of dimension 5x5x8 mm'^ cut along the crys- 
tallographic axes a, b and c (L = 8 mm) were used for 
these measurements. The samples were subjected to lin- 
ear temperature ramps of 0.2°C/min in the range — 20°C 
- -hl00°C. 

In a given temperature interval dT the linear thermal 
expansion coefficients are defined as 



1 dL, 



(4) 



(3) 



where Lt is the sample length along the direction i, with 
i = a,b, c. To be consistent with our assignment of the 
principal optical axes X, Y, Z in Section lVll we shall once 
and for all make the substitutions X = b,Y = a, Z = c 
in the following. 

For the measurement of the thermal expansion the 
light source of the absolute interferometric dilatometer 
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FIG. 3: Thermal expansion coefficients of LIS along the axes 
X = b,Y = a,Z = c from -20°C to +100°C. 



was a frequency stabilized 2-mW He-Ne laser. A plot of 
the experimental data is reported in Fig. 13 A linear fit 
of the data gives the following temperature dependence 
of the Unear thermal expansion coefficients, 



ax ^ 1.61 X 10' 
ay = 0.89 X 10" 
az = 0.66 X 10" 



+ 1.4 X 10"**r, 
ay = 0.89 X 10"^ + 0.7 x lO^^T, 



0.9 X 10"**r, 



(5) 
(6) 
(7) 



where T is the temperature (°C). The deviation of ex- 
perimental values from these fits are less than 10^^ K^^. 
At T = 20°C we have ax = 1.64 x lO-^K"!, ay = 
0.90 X 10-^ and az = 0.68 x IQ-^ K-\ respectively. 
Such values are of the same order of magnitude (or even 
less) than those for AGS {ai\c = 1.25 x 10"^ K"\ a±c = 
-1.32x10-5K-1) (ill and AGSe {a^^ = 1.68x10-5K-\ 
a^c = -0.78 X 10-5 K-i) jH. 

However, the main difference with those chalcopyrites 
is that LIS does not display their anomalous thermal 
expansion behaviour (opposite sign of a along the two 
orthogonal directions) that is a major source of thermo- 
mechanical stress in high power applications [s^ ■ In con- 
trast, LIS expands in the same direction as temperature 
increases and exhibits a rather weak temperature depen- 
dence. This favorable situation is related to the specific 
nature of the weak Li-S bond already evoked in the pre- 
vious subsection. 



C. Thermo-optic coefficients 

We shall again make use of the labels X = b,Y = 
a, Z = c denoting the principal optical axes to describe 
the change of the index of refraction versus temperature, 
for a light polarized along one of these principal axes 
and propagating along another principal direction. For a 
light wave polarized along direction pj and propagating 
along direction e^, we define the coefficient jji related to 
thermal changes in optical thickness (path) by 



1 din j Li) 
•i Li dT 



(8) 



where rij and L,; are respectively the refractive index and 
the length of the sample, with i, j = X, Y or Z and i ^ j- 



TABLE I: Principal normalized thermo-optic coefficients 
^ ^ ai + a2T of LIS [Eq. lEU] at four laser 

wavelengths. 



A (pim) 


A- (°c-^) 


lO'ai (°C-') 




0.4765 


Pz 


3.41 ± 0.06 
4.11 ± 0.07 
4.08 ± 0.02 


3.1 ± 0.3 

4.6 ± 0.9 

4.7 ± 0.4 


0.6328 


Px 
Py 
Pz 


2.28 ± 0.07 
2.76 ± 0.05 
2.69 ± 0.01 


1.5 ± 0.4 
2.1 ± 0.4 
2.4 ± 0.1 


1.0642 


pY 

Pz 


1.75 ± 0.07 
2.10 ± 0.11 
2.03 ± 0.08 


1.0 ± 0.06 
1.3 ± 0.5 
1.3 ± 0.1 


3.392 


Px 
Py 
Pz 


1.50 ± 0.07 
1.83 ± 0.14 
1.73 ± 0.03 


0.7 ± 0.17 
1.1 ± 0.5 
1.0 ± 0.02 



The normalized thermo-optic coefficients j3j = are 
derived from and © and given by 



Pj = 7ji 



(9) 



As evident from Eq.Q, the knowledge on the values of 
the thermal expansion is prerequisite for the interfero- 
metric determination of the thermo-optic coefficients at 
a given wavelength. For the orthorhombic symmetry of 
LIS it can be seen that each Pj may be obtained for two 
directions of light propagation. This has been performed 
at four laser wavelengths, and the resulting linear fits to 



Pj{\T) ^ ai{\) + a2{\)T 



(10) 



where T is measured in °C are summarized in Table 
The indicated uncertainties in the parameters ai and 
02 represent the deviation from the mean value de- 
termined from the two possible propagation directions 
Ci for a given orientation pj of light polarization. At 
A = 1064 nm and T = 20° C, one has dnx/dT = 3.725 x 
10-5 K-i, dny/dT = 4.545 x 10~^K-^, dnz/dT = 
4.467 X 10~5k~^. Compared with the thermo-optic co- 
efficients of AGS H^l, these values are about 5 times 
lower. Dispersion relations based on those wavelength- 
dependent data will be derived in Subsection IVI Bl 



IV. VIBRATIONAL PROPERTIES 

The vibrational properties have been investigated by 
means of Raman spectroscopy and infrared reflectivity in 
polarized light at room temperature. The group theory 
analysis shows that, at the zone centre for the /3-NaFe02 
structure, the 48 normal phonon modes are distributed 
among the various irreducible representations of the 



factor group as follows: 
F"'' = 12Ai - 



I2A2 + 12Bi + I2B2. 



(11) 



The three acoustic phonon modes have Ai, Bi and 
B2 symmetries; the 45 remaining (optical) modes are 
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Raman-active. Ai modes will be active when the experi- 
mental configuration selects a diagonal aa, bb or cc com- 
ponent of the Raman tensor, while the selection of ba, ca 
and cb components will respectively allow the observation 
of the optical phonons of A2, Bi and B2 symmetry. A2- 
type phonons being infrared inactive, the IR-reflectivity 
spectra should then contain 11 polar modes of each Ai, 
Bi and B2 symmetries, with dipole moment respectively 
parallel to c, a and b crystallographic axes. In theory, 
every optical phonon may then be observed by means 
of IR-reflectivity and/or Raman-scattering. Nonetheless, 
depending on various factors such as the more or less 
polar character of bonds, even an active mode may be 
unobservable. 

To obtain results on every phonon type, we used two 
LIS rose-tinge samples cut from the same ingot, with 
sizes 4x4x5 mm^ and parallel optically polished facets 
(4x4 mm^) normal to a and c axes. Infrared reflection 
spectra were recorded in the spectral range 18.3 - 600 
cm~^, and this whole range of data was used to adjust 
the parameters of the dielectric permittivity model shown 
in Ea. (|12|l . The experimental spectra, recorded at quasi- 
normal incidence (at an angle of incidence of 6°) with 
an electric field polarization parallel to a, b and c axes, 
are shown in Fig. ^ (symbols). As can be seen from 
these graphs, each reflectivity spectrum is rather complex 
and exhibits several more or less overlapping bands, and 
the high-energy band(s) in each polarization generally 
show some additional structure (slight inflexion in the 
band shape) that suggests, in addition to a strong mode, 
the presence of several other modes contributing to the 
same band. Actually, from the numerical analysis of the 
experimental results, it is found that for each direction 
of polarization, there are three polar modes in the high 
frequency band (350 - 420 cm~^) and medium frequency 
band (250-350 cm^^) and the remaining five polar modes, 
not all observable in Fig. 21 are in the low frequency 
region (< 250 cm~^). 

Polarized micro-Raman spectra were recorded using 
a Raman spectrometer equipped with a confocal micro- 
scope and a liquid-N2 cooled CCD detector in the region 
12 - 600 cm^^ in triple substractive configuration with an 
1800 gr/mm grating. To increase the quality and resolu- 
tion of the spectra, we chose to use the 647nm red line of 
an argon-krypton laser (P^lOO mW) which experiences 
lower absorption loss than a green laser ^J)- To avoid 
possible infiuence of surface defects, the confocal micro- 
scope focused the laser beam at a depth of 60^m under 
the surface. A scattering geometry will be described by 
the notation ei{piP2)e2 45] where identifies the prop- 
agation direction of the laser light, polarized along pi, e 
and p being a, 6 or c crystallographic axes; the subscripts 
i = 1 and 2 refer to incident and scattered beams respec- 
tively. Fig. 13 shows some of the experimental spectra. 
The Raman spectra are dominated by a strong mode of 
Ai symmetry at ^ 268 cm^^, also apparent in the IR 
reflectivity spectra for E//& and E//c. 

Due to the respective mass of the atoms, high (low) fre- 
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FIG. 4: Infrared reflection curves (symbols) recorded with 
polarized light with an electric field parallel to the crystallo- 
graphic axes a - Bi modes -, & - B2 modes - and c - Ai modes. 
The continuous lines represent in each case the best-fit reflec- 
tivity curves calculated with the model in Egs. 1121 - 1131 . 



quency modes obviously imply movements of Li-S (In-S) 
bonds. Phonons of medium frequency (250-350 cm~^), 
visible in both IR-reflectivity and Raman scattering, are 
attributed to In-S bonds ^3] while those of highest fre- 
quency (above 350 cm^^. Fig. 0)) that give rise to strong 
IR bands but no observable Raman line are due to weaker 
Li-S bonds Js^l . Polar modes TO«LO splitting observed 
in the infrared measurements for any of the explored po- 
larizations is rather important in the high frequency re- 
gion (strong IR-bands) which attests for a high polar 
(ionicity) character of Li-S bonds in this crystal. The 
splitting decreases with the phonon frequency to become 
so small for low frequency modes that no corresponding 
band appears in the reflectivity spectrum which indicates 
a mostly covalent character of the In-S bonds. However, 
even if a phonon cannot be observed by means of IR- 
reflectivity, its Raman activity may be high, which is the 
case for modes lying below 100 cm^^. Actually intense 
Raman lines appear in the low and mid-frequency region, 
when almost no line can be detected above roughly 350 
cm~^. In order to flt the experimental IR-reflectivity 
spectra, and get the best-flt reflectivity curves (lines in 
Figure^ in a good agreement with experimental points) 
we used the so-called four-parameter model, convenient 
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FIG. 5: Polarized Raman spectra (shifted for reasons of clar- 
ity) recorded at room temperature. The experimental config- 
urations b{aa)b, c{bh)c, b{cc)b correspond to Ai modes, c{ba)c 
to A2 modes, b{ca)b to Bi modes and a{cb)a to B2 modes. 



to represent even asymmetrical and wide bands |4^| . This 
model expands the factorized form of the dielectric per- 
mittivity function as 



N 



e{ijj) = e - is = £00 



- -I- ijj^TO 



1 "j,TO 



(12) 



to fit a reflectivity curve described by the Fresnel formula 



(13) 



In the /3-NaFe02 structure, the model implies for a given 
polarization N damped oscillators, the number N be- 
ing at most equal to 11 (number of theoretically observ- 
able modes, predicted in this structure by group the- 
ory analysis). The j-th damped oscillator represents the 
j-th observed infrared active couple of modes, charac- 
terized by four adjustable parameters, namely the lon- 
gitudinal rJj.LO and transverse flj.TO angular frequen- 
cies, and the corresponding damping frequencies 7j,LO 
and 7j,TO- The high-frequency limit values of the di- 
electric constant is also adjusted to yield £oo(Ai)=3.99, 
eoo(Bi)=4.09, eo,(B2)=3.97. 



Tablelnllists the wavenumbers 



ing rates 7^ = jj/2'Kc (expressed in cm ) of the observed 



rij/27rc and damp- 



optical phonons deduced from IR-reflectivity and/or Ra- 
man scattering experiments, assigned to the various ir- 
reducible representations of the C^y factor group. These 
data allow to compute the linear dielectric properties of 
LIS in the far-IR range, using Ea. ((T^ . Typically, TO and 
LO frequencies deduced from IR-reflectivity experiments 
are determined to better than 1 cm~^. The precision of a 
Raman-active phonon frequency is also generally better 
than 1 cm~^, except for broad and low- intensity lines for 
which it decreases to 2 cm~^. The transverse or longi- 
tudinal character of the Raman-active modes for Ai, Bi 
and B2 symmetries has been determined with the help of 
the IR reflectivity results where this TO or LO character 
is immediate. Non-observed IR-active modes indicate a 
small TOwLO splitting and generally did not allow an 
accurate determination of TO and LO frequencies from 
Raman spectra, but only of a mean value labelled by 
TO«LO. A few Raman lines have been attributed to an 
interference by an intense line active in another geome- 
try. These Raman results complement those fromprevi- 
ous studies performed with unpolarized l26l |23 | and 
polarized light B3 | . 

The accurate determination of the frequencies in Ta- 
ble ^ allows to estimate the interatomic bonding force 
constants in LIS, using the theoretical considerations de- 
veloped by Neumann about the correlation between the 
vibrational modes of compounds with sphalerite, chal- 
copyrite and /3-NaFe02 structures According to 

this model the vibrational spectrum of a crystal with /3- 
NaFe02 structure contains two sphalerite-like modes the 
TO-frequencies of which are nearly independent of the 
polarization direction. Further, one of these sphalerite- 
like modes always corresponds to the highest Dto fre- 
quency found in the IR reflectivity spectra. If the in- 
teraction between the atoms is restricted to next-nearest 
neighbors which is justified by the high bond ionicities 
in Li -containing ternary compounds, the wavenumbers 
of the transverse optical modes of these sphalerite-like 
modes in LiB^^^C^^ are given by 



1 



TO,k 



■{ak ~ Sak) 



1 

TOfc 



1 

mc 



(14) 



with the superscript fc = in the first parenthesis denot- 
ing Li-C or B-C cation-anion bonds (in the present case, 
C —S and B =In). Here c is the speed of light, 
the mechanical bond-stretching force constant, ruk the 
mass of the cation and mc the mass on the chalcogen 
anion (C^S^^). The quantity Sak accounts for the fact 
that in the case of IR active modes the local atomic dis- 
placements give rise to a non-vanishing net dipole mo- 
ment per unit cell and thus, to additional dipole-dipole 
interactions. It is given by Sat = V3(efc)2/(16eos^), 
where Sk is the average bond length {sLi-s — 0.2438 
nm, sin-s = 0.2452 nm), Ck = 0.28e/fe(Zfe + Zc) is 
an effective charge depending on the bond ionicity fk 
and the valences of the cation and anion participating to 
the bond. In Ref. [l^, from an unpolarized IR reflectiv- 
ity spectrum of polycrystalline LIS restricted to the re- 
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TABLE II: Zone-center phonon wavenumbers (S'to/S'lo) and damping (7T0/7L0), in cm~^, deduced from IR reflectivity 
and Raman scattering experiments at room temperature for the various symmetry-adapted polarization conflgurations. The 
transverse (TO) and/or longitudinal (LO) character of the observed Raman-active phonons have been determined with the 
help of the IR reflectivity results (the notation TO~LO means that the splitting is too weak to be unambiguously identified). 



Phonon 


IK reflectivity 




Raman 


symmetry 


Wavenumbers Uj 


Dampings 7, 


scattering 


Ai 
TO/LO 
(cm-i) 


t 

128.6/138.2 
266.0/267.1 
274.2/279.2 
274.2/279.2 
313 6/313 8 
323.7/333.7 
341.7/356.0 
357.0/373.0 
373.1/399.0 


3.0/4.4 
8.5/8.7 
9.2/8.7 
9.2/8.7 
4.9/5.0 
6.2/11.3 

14.3/13.9 
14.5/14.7 
14.3/16.5 


56.4 (TOf«LO) 
88.1 (TOr^LO) 
112.7/113.3 
129.1/135.2 
268/269 
273/- 
273/- 
314 fTO^iLO") 
323.4/333 
341/- 
357/- 
-/399 


A2 
(cm ) 


IR inactive 




59.4 
68.7 
117.8 
125.9 
157.4 
253.5 
297.4 
319.4 
332.5 


Bi 
TO/LO 

(cm-i) 


125.8^133.4 
164.9/166.0 
264.9/266.3 
303.3/303.7 
318.3/334.6 
343.6/347.8 
356.4/375.9 


2.6^/3.0 
8.4/9.0 
9.3/10.3 
9.4/9.2 
5.0/7.2 
11.2/10.4 
10.4/22.8 


61.1 (T0«L0) 
93.5 (TOf«LO) 

125.8/133.3 
165 (TO^LO) 
267(TO^LO) 

3187/ - 

356^5/ - 



376.7/395.8 21.9/23.1 -/- 
400.4/408.6 29.7/12.4 -[-_ 









59.3(TOf»LO) 




96.1/96.4 


3.3/3.1 


96 (TOf^LO) 




104.7/105.7 


4.8/5.2 


104.7/105.5 


B2 

TO/LO 
(cm-i) 


132.0/132.5 


10.2/10.1 


133 (TO^LO) 


142.7/144.9 


6.1/6.9 


143.2 (TOsiLO) 


276.5/285.1 


8.0/8.4 


276.2/285 


324.6/329.8 


5.9/7.3 


325.2 (TO^LO) 




333.4/336.0 


7.2/7.5 


334 (TO^LO) 




348.0/358.4 


10.3/8.9 


348/- 




362.2/397.7 


8.6/11.0 


362/- 




398.6/410.4 


11.1/11.9 


-/- 



gion 200 - 400 cm ^ , the highest-frequency TO /LO mode 
identified corresponds to 356/406 cm^^ and the second 
TO/LO sphalerite- like mode to 323/334 cm"^ Actually 
from the inspection of Table (nj the highest-frequency 
sphalerite-like modes with nearly identical frequencies in 
Bi and B2 symmetries are respectively the 400.4/408.6 
cm~^ and 398.6/410.4 cm~^ (this mode is not observed in 
Ai). The second sphalerite-like mode is recognized from 
Tableinias the 323.7/333.7 cm-i (Ai), 318.3/334.6 cm-^ 
(Bi), 324.6/329.8 cm~i (B2), close to the value identified 
in Ref. "3^. The Li-S bond force constant deduced from 
Eq.(in|) by Sobotta et al an = 24.7 N/m, is then 
underestimated by (400.4/356)^ = 24% while the In-S 
bond force constant, deduced from the second sphalerite- 
like mode remains unchanged to a/„ = 53.4 N/m. But 



even with this slight correction, the Li-S bond remains 
still about twice weaker than the In-S bond. In con- 
trast the strength of the Ag-C^' bond in chalcopyrite 
compounds is more comparable to that of the B^^^-C^^ 
bonds jS^ , underlining the specific nature of the lithium- 
chalcogen bond in crystals with tetrahedral coordination 
of the atoms. As a consequence, the increased lattice 
phonon energy of LIS, as compared with AGS or AGSe, 
favors heat dissipation, hence a higher thermal conduc- 
tivity is expected. Indeed, an unpublished measurement 
of LIS thermal conductivity yielded Ka = 6.0 W/(m-K), 
Kb = 6.2 W/(m-K) and = 7.6 W/(m-K) 13, i.e. 5 
times larger than for AGS and of the same magnitude as 
for LiNbOg. 

A last important feature from the vibrational spectra is 
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that the phonon spectrum for LIS is located at wavenum- 
bers below 410 cm"^. The IR wavelength cut-off (8.1 
/im) in Fig. 13 quite exactly corresponds to three times 
this highest energy phonon. The long wavelength shape 
of the transmission curves - and limit of the transparency 
range - seems then to be simply due to multiphonon pro- 
cesses. It is worth mentioning that, although the vibra- 
tional investigations have used rose tinge samples, the 
phonon spectrum should not depend on coloration since 
all LIS samples display the same IR cutoff edge, with or 
without annealing. 

V. PIEZO-ELECTRIC AND ELECTRO-OPTIC 
COEFFICIENTS 

A. Apparent electro-optic coefficients 

The measurement of electro-optic coefficients involves 
the determination of changes in optical thickness when 
applying an electric field to the material. LIS belongs 
to the symmetry class mm2 for which the electro-optic 
tensor exhibits 5 nonzero components, namely ri3, r23, 
?'33, ?'42 and rsi [s^. In this subsection we will describe 
the electro-optic effect in the principal optic frame where 
XYZ are the axes of the index ellipsoid which as we shall 
see in the next section coincide with the crystallographic 
bac axes. We will address only the r^a's coefficients that 
can be obtained by applying a field E^, along the Z-axis 
which does not lead to a rotation of the ellipsoid. Opti- 
cal interferometric methods may be advantageously used 
in this case for their determination. An applied electric 
field E, with respective components {Ei, E2, E3) along 
the principal axes {X, Y, Z) , will modify the index ellip- 
soid according to 

+ inf + r33E3)Z^ + 2ri2E2YZ + 2rr,^EiXZ = 1. (15) 

The new indices riilE^) (i = 1, 2, 3 = X, Y, Z) are given 
by 

n^{E3) = - ]^nyaE3 (16) 

where 5{ni) = —^nfri^E^ represents the variation of the 
refractive index n^. 

Due to the piezo-electric effect the field E3 will also 
generate a variation 5{Lj) of the crystal thickness Lj 
along a given direction ej such that 

S{L,) = L,d3,E3 (17) 

where d^j {j — 1,2,3) are components of the piezo- 
electric tensor ^^j. For a light wave propagating along 
direction Cj and polarized along direction pi of the crys- 
tal the variation S{niLj) of the optical path UiLj induced 
by an applied electric field E3 will be 

5{niLj) = UiLj (dsj - ^nfrisj E3. (18) 



TABLE III: Apparent electro-optic coefficients of LIS; 3(X) 
and 3(Y) represent the values of rJjf obtained for a light 
propagation along X and along Y, respectively. 



i 


r^ff (pm/V) 


1 


2.18±0.04 


2 


2.82±0.33 


3(X) 


0.93±0.03 


3(Y) 


0.11±0.02 



We define the effective electro-optic coefficients as 

r'i = -^^{d:i,-\nln^ (19) 

with j = 1,2, which represent in fact the (apparent) 
electro-optic coefficients obtained in most experiments 
where the piezo-electric contribution is ignored. The 
measurements of the r^^'s were performed by using the 
technique of Phase-modulated Thermal Scanning Inter- 
ferometry (PTSI) [sj. In this method, as for the mea- 
surements of thermo-optic coefficients, the sample acts as 
a Fabry-Perot thermal scanning interferometer in which 
phase modulation is achieved by applying an AC field E3 
along the Z axis. Various configurations of light polar- 
ization {Ex,Ey,Ez) and wave vector k {kx,kY) were 
used to determine the apparent electro-optic coefficients 
of LIS. 

The two LIS samples of aperture 5x5 mm^ used for the 
thermo-optic data measurements with lengths L = 8 mm 
along the X and Y propagation axes were supplied with 
gold evaporated electrodes on their Z facets. The low fre- 
quency (480 Hz) applied voltage was 200 V peak-to-peak 
and linear temperature ramps of 0.2°C/min were used 
from 20°C to 55°C. Measurements of apparent electro- 
optic coefficients were performed at 1.064 /im by using a 
lOOmW TEMqo single-frequency Nd:YAG laser. In the 
scanned temperature interval a constant value has been 
observed for each of them; they are summarized in Ta- 
ble mil These data feature a quite weak electro-optic 
effect in LIS and the large discrepancy observed between 
^i3(x) ^^'^ '''i3{Y) suggests a strong contribution of the 
piezo-electric effect to the apparent electro-optic one. 

B. Piezo-electric and pure electro-optic coefficients 

The non- vanishing piezo-electric coefficients ^31, ^32 
and ^33 were measured by using the modified Mach- 
Zehnder arrangement of the absolute interferometric 
dilatometer in a manner similar to the one employed for 
dilatation measurements. Phase modulation of the inter- 
ference pattern is performed by applying to the samples 
the same voltage and same thermal scanning procedure 
as described previously. The algebraic sign of the co- 
efficients was attributed according to the ANSI/IEEE 
Standard 176-1987 on piezo-electric crystals 36] and 
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we have found in pm/V units and in the XYZ frame: 
= -5.6 ±0.1 , J32 = -2.7 ±0.2 and (I33 = +7.5 ±0.2. 
From Equation (|19|l we determine the pure electro- 
optic coefBcients ri3 = 1.0 pm/V, r23 — 0.4 pm/V and 
r33 = —1.3 pm/V which reveals a rather small true 
electro-optic effect compared to the piezo-electric one. 
The relative uncertainty in these values is close to 15%. 
These measurements confirm that, as AGS or AGSe, 
LIS is a poor candidate for electro-optic or piezo-electric 
based devices. A comparison of our results with ear- 
lier data of Negran et al [S^l cannot be meaningful since 
the ri3, r23 and r33 that were measured there did not 
account for the piezo-electric contribution and were not 
determined independently, but as a linear combination. 



requires at first an accurate determination of the piezo- 
electric contribution. 



VI. LINEAR OPTICAL AND THERMO-OPTIC 
DISPERSIONS 

We now address, in light of the vibrational spectra re- 
sults, the most important data for a nonlinear material: 
linear optical and thermo-optic dispersion relations that 
allow to predict accurately the phase-matching directions 
and temperature tunability. 



C. Comparison of the electro-optic and 
pyroelectric coefficients 

From the anharmonic oscillator model developed by 
Garrett [S^, Soref [E^l yields a semi-empirical rela- 
tionship between the pyroelectric coefficient at constant 
stress {p'^ — and the second-order electronic 

susceptibility x'^^ describing the linear (true) electro- 
optic effect: 

(20) 

The pyroelectric coefficient p"' gives the variation of the 
spontaneous polarization APs induced by a temperature 
variation AT. The parameters and e are the mass and 
charge of the electron, uj^ and oji are respectively the elec- 
tronic and ionic resonance frequencies, Cy is the molar 
specific heat at constant volume (which differs from the 
specific heat at constant pressure Cp measured in Subsec- 
tion llll Al bv a negligible quantity) and Na the Avogadro 
number. The constant terms A, B, C and D are the coef- 
ficients related to the third order expansion of the crystal 
potential energy in the electronic and ionic coordinate 
system qe and 54]. The ratio (SA + B)/{C + 3D) 
indicates the extent to which the predominant ionic con- 
tribution to the anharmonic potential at w uji differs 
from the predominant electronic one at uj >> uii. This 
ratio is expected to be constant whatever the material, 
suggesting a linear relationship between the electro-optic 
and pyroelectric coefficients. This hypothesis was verified 
for a few known materials, mostly of the oxide type ^5^. 

By using Cy = Cp and X333 = 7T,|r33 from the previ- 
ous subsection, the pyroelectric coefRcient measured pre- 
viously "551 as a function of temperature and its value 
extrapolated at 300 K, p'" = 6 x 10^1° C- cm-^K-\ 
we observe that Soref's hypothesis is confirmed for the 
ternary semiconductor LIS, as it was already shown for 
the binary one CdS [H^l- We note also that a rather 
good estimate of the pyroelectric effect may be deduced 
from electro-optic measurements and vice versa, but this 



A. Sellmeier equations at room-temperature 

The only available dispersion data of LIS, measured 
using the minimum deviation technique, were compiled 
by Boyd et al in 1973 and fitted to one-pole (UV) Sell- 
meier equations ^5^ reproduced also in Refs. |2lll5^l58j |. 
Although the accuracy of their index determination was 
not reported, from subsequent nonlinear phase-matched 
conversion experiments it appeared that the Sellmeier 
equations based on their data are accurate to predict the 
birefringence (via angle acceptance bandwith measure- 
ments [13): but fail to predict accurately the absolute 
phase-matching directions. Discrepancies of several de- 
grees from the expected phase-matching angles were re- 
ported [23 ■ Such deviations cannot be attributed to the 
X-ray sample orientation uncertainties (±0.5°). Further- 
more, we have checked that the phase-matching angles, 
and hence the linear refractive indices, are independent 
of crystal coloration or exact stoichiometry. 

Following the standard convention for assigning the 
principal axes {X, Y, Z) to the crystallographic axes 
(a, 6, c), e.g. nx < ny < nz, one has X ^ b, Y ^ a and 
Z <-!■ c for LIS since nb < na < tie. A phase-matching 
direction will be labelled by {9, cp), where 6 is the angle of 
the wavevector respective to the polar Z axis and (p is the 
azimuthal angle between its projection in the X — Y plane 
and the X axis. To refit the Sellmeier equations given in 
Ref. IH^, we have used the experimental type-II (eoe, 
oeo) SHG angles measured in the fundamental wave- 
length range A e [2.5 — 6] /xm 113,0], supplemented with 
new recent data derived from type-II (eoe) difference- 
frequency generation (DFG) in the X-Y plane (see sub- 
section IIX B|) . where in the wavelength-tuning mode, 
mid-IR down-conversion in the range Ai G [6.6 — 7] /im 
was achieved at normal incidence. Table Hvl compiles the 
experimental phase-matching data used to derive the fol- 
lowing room-temperature (T = 20°C) two-pole Sellmeier 
equations, valid over the full transparency window, using 
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a nonlinear least-square fit of the data, 



njf = 6.686059 - 



7.095493 



7.256327- 





0.1385833 


2047.46509 


A2 


- 0.05910334 ^ 


A2 - 897.7476 




0.1422326 


2511.08936 


A2 


- 0.06614640 ^ 


A2 - 988.2024 




0.15072 


2626.10840 


A2 


- 0.06823652 


A2 - 983.0503 



,(21) 
(22) 
(23) 



In Eas. (|21|l - H23(l . the wavelengths are expressed in 
microns. The three IR wavelength poles Xir{X) = 
29.96 /im, XiB.{Y) = 31.44 ^m, Xir{Z) = 31.35 /im 
(i.e. D{X) = 333.7 cm-i, D{Y) = 318.1 cni"\ D{Z) = 
318.9 cm^^) correspond to the medium-frequency phonon 
band-centers of symmetry B2, Bi and Ai respectively, at- 
tributed to the In-S vibrations (see Table llT|l . We note 
that another two-pole Sellmeier set of equations that re- 
cently appeared |53| is based on the old index data of 
Boyd et al. |^ while the mid-IR poles were not associated 
with phonon frequencies. Fig. |S1 displays the principal in- 
dices of refraction deduced from Eas. lf^ - lf^ . No index 
crossing (no(A) = ne(A)) is found in the blue part of the 
dispersion data, in contrast with the uniaxial AGS [43l |. 
The two optic axes, lying in the X-Z plane, determine 
the propagation directions where the index of refraction 
is independent of the polarization. The angle between 
the two optic axes is 2Vz where Vz is the angle of each 
axis with the polar Z-axis. Vz is wavelength-dependent 
and can be calculated from Q 



sinV^ 



nz{n 



AY"' 



(24) 



Fig. [3 plots Vz as a function of the wavelength, using 
Eas. H21|l -H23 |l (solid line) and using the Sellmeier equa- 
tions of Ref [53 (dashed line). Let us note the different 
behavior for A < 1 /im (absence of a retracing behavior). 
Note that although the principal indices on which the 
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FIG. 6: Principal refractive indices of LiInS2 as computed 
from Eqs.l|2IJ-lIlSl- The inset plots show the UV and mid-IR 
portions on an expanded scale and also a comparison with a 
calculation based on the Sellmeier expansions from Ref. [H^ 
(dashed lines). 
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FIG. 7; Angle Vz between the optic axes and the Z-axis 
in LIS calculated with the Sellmeier expansion of Ref. [H^ 
(dashed line), and the new Sellmeier expansion (Eas. (l2H - 
^) (solid line). 



old dispersion relations are based were measured from 
the minimum deviation technique in 1973 'oj, the ab- 
sorption coefhcient of the prism samples used then was 
quite different (much larger than in Fig. [5] up to about 
2.5 /im) which obviously is the reason for the difference 
observed in Fig. [7| Since Vz > 45° LIS is a negative 
biaxial crystal. At longer wavelengths it can be regarded 
as quasi-uniaxial, which is a consequence of the fact that 
the two indices ny and nz are very close. Therefore no 
phase-matching can be expected for propagation direc- 
tions in the vicinity of the X-principal axis. 

The agreement between the experimental phase- 
matching data used to derive Eas. H21|) - (|23|l and the cal- 
culated ones is illustrated in Fig|Sl (solid lines) and Ta- 
ble IIVI The accuracy of phase-matching angle predic- 
tion is ±1.5°. The eight last SHG points in Fig. |S1 were 
discarded from the nonlinear fit for obvious lack of relia- 
bility. Only a two-pole form leads to the smallest resid- 
ual discrepancies between calculated and experimental 
data, while the former one-pole Sellmeier equations S9| 
are unable to predict SHG near the cut-off wavelength 
of ~ 6 /tm. Let us note that there is a discrepancy of 
~ -|-6° between the experimental phase-matching angles 
for SHG in the Y-Z plane and the ones calculated with 
the old Sellmeier equation, which means that the Y-Z 
birefringence given by Eqs. (|21|I - H23|I is more accurate. 
The previous one-pole Sellmeier equations also predict 
DFG wavelengths ~ -1-40 nm longer than experimentally 
observed. Under the usual designation Ai > A2 > A3 the 
FWHM DFG angular acceptance (FWHM) measured at 
fixed wavelengths was = 0.19° for a crystal length of 
10 mm, i.e. much smaller than the uncertainty on crys- 
tal orientation (±0.5°). Hence the uncertainty (±1 nm at 
maximum) of wavelength measurement in Table ITVl can- 
not explain the ~ 40 nm shift observed with the previous 
Sellmeier equations, confirming the increased accuracy in 
dispersion provided by Eas. H21|l - H23|) . 
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TABLE IV: Phase-matching conditions for type-II SHG and DFG in LiInS2. 



Wavelength (nm)" Phase-Matching Angles {0,<^Y 





Ai 


A2 


A3 iVieasured 


Gaicuiated 


SHG 


2366 


2366 


1183 


90,82.1 




90, 83.908 




(XY: eoe) 


2469 


2469 


1234.5 


90, 73.1 




90, 72.755 




2527 


2527 


1263.5 


; 90, 69.8 




90, 69.235; 






2583 


2583 


1291.5 


'90, 67.4 




90, 66.507, 






2900 


2900 


1450 


90,57.9 




90, 56.884 






3400 


3400 


1700 


; 90, 50.7 




90, 50.543; 






3700 


3700 


1850 


'90, 48.3 




90,49.392; 






3900 


3900 


1950 


90, 49.0 




90,49.381 






4450 


4450 


2225 


; 90, 51.6 




90,51.772; 






4950 


4950 


2475 


'90,56.3 




90, 56.543) 






5350 


5350 


2675 


; 90, 62.8 


) 


;90, 62.227; 






5550 


5550 


2775 


'90, 66.0 




'90,65.916) 






5750 


5750 


2875 


90,69.1 




90, 70.508 






5900 


5900 


2950 


;90, 71.4 




;90, 75.006, 




SHG 


2542.7 


2542.7 


1371.35 


'35.4, 90 


) 


'36.452, 90 




(YZ: eoo) 


2552.7 


2552.7 


1276.35 


;34.0, 90 




;34.684, 90 




2570.4 


2570.4 


1285.2 


'31.0,90 




'31.508, 90 






2587 


2587 


1293.5 


28.7, 90 




28.435, 90 






2602.3 


2602.3 


loUl.io 


oc; on 




<op: /1 01 on 






/DUD. ( 


2606. / 


1303.35 


'25.1,90 




'24.592, 90 






2631.4 


2631.4 


1315.7 


19.6,90 




'19.211,90 




DFG 


6585.14 


872.240 


770.220 


(90, 42) 




'90,42.170) 




(XY: eoe) 


6623.95 


865.801 


765.716 


(90, 42) 




;90, 42.145) 




6719.87 


850.443 


754.905 


(90, 42) 




'90,42.076) 






6762.44 


843.761 


750.162 


(90, 42) 




90,42.046 






6832.63 


833.187 


742.629 


(90, 42) 




;90, 41.992) 






6899.92 


822.391 


734.810 


(90, 42) 




'90,41.967) 






6989.56 


810.827 


726.544 


(90, 42) 




;90, 41.863) 





"l/Ai -I- I/A2 = I/A3 with Ai > A2 > A3, 
'in degree. 
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FIG. 8: Phase-matching angles for type-II (eoe or oee) SHG 
in the X-Y principal plane (main frame) and type-II (oeo or 
eoo) SHG in the Y-Z plane (inset frame) of LIS. The funda- 
mental source for A < 2.65 /^m (A > 2.75 /xm, respectively) 
is the idler wave of a Nd:YAG-pumped nanosecond LiNb03 
parametric oscillator (the frequency-doubled radiation from a 
Free- Electron Laser). The symbols refer to experiments while 
the solid lines are calculated from Eas. ll2H - l231 . The dashed 
hues are computed from the Sellmeier equations of Refs. [H^. 



B. Thermo-optic dispersion relations 

The thermo-optic coefficients /3j [Eq. (|10() ] are func- 
tions of both temperature T (in °C) and wavelength A 
(in \im). In analogy with the index dispersion equations 
at room-temperature (see previous subsection) the pa- 
rameters ai;(A) (k = 1,2) in Tableware fitted according 
to the following two-pole functional form, 



afe(A) = Co,fe + 



2.k 



A2 



^01, fe 



A2 



^02, fe 



(25) 



Table fVl compiles the fitting parameters Co,fc, Ci.fe, C2,fc, 
Aoi,fc and Ao2,fc- According to the definition Pj = 
the index of refraction at a given temperature T will be 
given by 



n{X,To) 



exp 



ai(A)(r-To) + a2(A) 



(r-ro)2i 



1 + ai(A)(r - To) + a2(A)^^— ^ (26) 



where Tq = 20°C and n(A,ro) is given by Eqs.^-^. 

These thermo-optic dispersion relations will have to 
be refined by experimental temperature phase-matching 
data in the future. Relations H25|) with the constant val- 
ues in Table would provide a suitable basis for such 
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TABLE V: Fitting parameters describing the thermal and wavelength dependence of the principal thermo-optic coefficients 
of LIS, according to Eas.(??)-fEii. 



afc(A) 


Pi 


G) 


Ci 




C2 


Aoi 


Ao2 




X 


2.899680 xlO"'^ 


0.2033761x10" 


-5 


297.0535x10"^ 


0.3475601 


14.86755 


ai 


Y 


2.478267x10"'^ 


0.2612463x10" 


-5 


122.3207x10"^ 


0.3391212 


13.91854 




Z 


3.315830x10"^ 


0.2555208x10" 


-5 


364.9813x10"^ 


0.3451539 


15.44113 




X 


1.987083x10"* 


0.1860292x10" 


-8 


250.0777x10"* 


0.3886392 


14.25991 




Y 


0.1523235x10"* 


0.2601335x10" 


-8 


-191.8712x10"* 


0.3969192 


14.79803 




Z 


-06564806x10"* 


0.4291068x10" 


-8 


-343.1015x10"* 


0.3435849 


14.94811 



a refit procedure. They predict a FWHM temperature 
acceptance of ATL — 248°C-cm (L: crystal length) at 
room temperature for the type-II (eoe) SHG at the fun- 
damental A — 2.5 /im in the X-Y plane. Hence LIS is 
expected to be very stable against self-induced thermal 
effects. As a counterpart, its temperature tunability is 
expected to be somewhat limited. However for the DFG 
of near-IR sources (A3 = 754.905 nm, A2 = 850.443 nm, 
see Table HVj) . the FWHM temperature acceptance is pre- 
dicted to be only ATL = 46°C-cm. 



VII. PHASE-MATCHING INVESTIGATIONS 

A. Effective nonlinearity and Hobden classification 

The phase-matching loci and the effective nonlinear- 
ity for LIS will be analyzed in the XYZ-principal op- 
tic axis frame, but the tensor elements of the non- 
linear susceptibility are defined traditionally, in accor- 
dance with ANSI/IEEE Std 176-1987 (3^, in the abc- 
crystallographic frame where c is the polar two-fold axis 
as in Ref. 9|- The chosen convention nx < ny < nz 
in order to have the two optic axes in the X-Z principal 
plane affects therefore the expressions for the effective 
second order nonlinearity doff. 

The general form of the contracted d^^^ tensor for the 
orthorhombic class mm2 reads as 



/ di5 
d(2) = ^24 

Vrfsi ^32 ^33 



(27) 



For an arbitrary propagation direction assuming 
coUinear interaction and neglecting the spatial walk-off 
effect analytical expressions for d^s can be derived in 
our conventions from the general formulae presented by 



Dmitriev et al l6C 



d^a^ = 2 di5 sin 9 cos 5 (cos sin ip sin 5 — cos tp cos 5) 



-'off 



X (cos 9 sin ip cos 5 + cos f sin 5) 

+ 2 d2i sin 9 cos 5 (cos 9 cos ip cos 5 — sin p sin 5) 

X (cos 9 cos If sin 5 + sin ip cos 5) 

+ dai sin 9 sin 5 (cos 9 sin ip cos 5 + cos ip sin 5^ 

+ d^2 sin 9 sin 5 (cos 9 cos p cos 5 — sin p sin 5)^ 



das sin"^ 9 cos^ 5 sin 5 



(28) 



and 



d^''^ - d 

"eff — "off 



= — di5 [sin 9 cos 6 (cos 9 sin pi sin 6 — cos p cos (5)^ 
-I- sin 9 sin 6 (cos 9 sin p) sin S — cos p cos S) 
X (cos 9 sin p) cos 5 + cos p sin 5)] 

— d24 [sin 9 cos S (cos 9 sin p sin S + sin p cos 5)^ 
+ sin 9 sin S (cos 9 cos p cos (5 — sin p sin (5) 
X (cos 9 cos (/? sin (5 -|- sin cos S)] 

— dsi sin 9 sin (5 (cos sin p sin (5 — cos p cos (5) 
X (cos 9 sin cos 6 + cos (/3 sin 5) 

— d^2 sin sin 5 (cos cos p cos (5 — sin p sin (5) 
X (cos 9 cos sin (5 -f sin (/s cos (5) 

— sin"^ 9 sin^ 5 cos (5 

(29) 



for type-I and type-II phase-matching, respectively, 
where the superscripts " s" and " /" stand for the " slow" 
and " fast" eigenmodes of polarization and their sequence 
follows the usual convention A1A2A3 with Ai > A2 > A3. 
The angle 8 whose introduction simplifies the expressions 
and which is determined from 



tan 25 



cos 9 sin 2p 



cot Vz sin 9 + sin p — cos-^ 9 cos^ p 



(30) 



(0 < 2(5 < tt) defines the polarization directions of the 
slow and fast waves which are orthogonal to each other. 
It is the angle between the polarization direction of the 
slow wave and the plane of propagation containing the 
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principal Z-axis. An approach for the generahzation of 
the doff expressions in the case of colhnear interactions 
taki ng in to account the spatial walk-off can be found in 
Ref. |6ll. 

In the principal planes the expressions for dcs l( ^ - (f^ 
are reduced to: 

jeoe ^oee 

"eff — "cff 

= -(^24 sin^ if + di5 cos^ tp) (X-Y plane) (31) 

d°^° ^ d^^" = -d24 sine (Y-Z plane) (32) 

= +d3i sine* (X-Z plane, 6* < Fz) (33) 

= 4^^° = -di5 sine (X-Z plane, 6* > Vz) (34) 

with superscripts " o" and " e" denoting the ordinary and 
extraordinary beams. LIS behaves as an optically neg- 
ative uniaxial crystal in the X-Y and X-Z (for e < Vz) 
planes and as an optically positive uniaxial crystal in the 
Y-Z and X-Z (for 6* > Vz) planes. Assuming the Klein- 
man symmetry condition to hold then di^ = c?3i and 

1^24 = ^32- 

The phase-matching loci for colhnear SHG can be cat- 
egorized using the classification of Hobden [g^I • Our cal- 
culations based on the relatively simple transcendental 
equations derived by Yao et al |63| , predict that no phase- 
matched SHG is possible below 1617 nm. LIS enters Hob- 
den class 13 at 1617 nm, class 11 at 1783 nm, class 10 
at 2353 nm, and class 9 at 2675 nm, then it goes back 
to class 10 at 5493 nm, to class 11 at 6111 nm, to class 
13 at 8224 nm, and finally no SHG phase-matching is 
possible again above 8710 nm. These transitional funda- 
mental wavelengths correspond to propagation directions 
along the Y - or along the Z - principal axes (noncriti- 
cal phase-matching) either for the ssf or the fsf {sff) 
polarization configurations. 

Figs. |HJi,b illustrate the different classes for SHG in 
LIS at several representative fundamental wavelengths. 
The surface of the unit sphere is projected onto the X-Z 
plane of the crystal. The direction of the wave vectors 
of the interacting waves for phase-matching as given by 
their interception with the surface of the unit sphere is 
plotted. The stereographic projection of the first octant 
is presented only but the loci in the other octants can be 
obtained by mirror reflections across the principal planes. 
We do not present plots of the angle S which can be easily 
calculated from Ea. (|Sn|l . Its dispersive properties (deter- 
mined by the angle Vz) lead to slightly different values at 
the three interacting wavelengths but normally the effect 
on the conversion efficiency can be neglected. Qualitative 
arguments relating the spatial walk-off (magnitude and 
direction) to the topology depicted in Fig.^can be found 
in Ref. |33- The calculation of the spatial walk-off for 
an arbitrary propagation direction is beyond the scope of 
this paper but several works provide suitable approaches 
to this aim [g^, EE Ez| • Analytical approaches for 
the estimation of acceptance parameters based on the 
small-signal approximation in the general case of biax- 
ial crystals can be found in Ref. |63, E^ • We note here 
only, that since the angle Vz for LIS is not far from 90° 




fS » 45 » 75 n 




IS W 43 « 75 n 



FIG. 9; (a)-(b): Stereographic projections of the SHG in 
the first octant of LIS calculated for wavelengths represen- 
tative of the Hobden classes. Type-I (ssf) interaction (solid 
lines) and type-II {fsf, sff) interactions (dashed lines). The 
wavelength annotations refer to the fundamental ones. 

the dispersive properties and in particular the spectral 
acceptance are not expected to be substantially modi- 
fied by propagation outside the principal planes. Fig. |H| 
indicates that double solutions for the angles exist (e.g. 
for Hobden's class 9) which is equivalent to the existence 
of points outside the principal planes where dtp/dO = 
holds. This means that noncritical phase-matching in 
one direction can also occur outside the principal planes 
(e.g. at = 26.7°, 6* = 74.35° for ssf type SHG at 4^m). 



B. SHG in principal planes 

The SHG phase-matching directions in the principal 
planes are shown in the lower part of Fig. 1101 The two 
first panels (X-Y) and (Y-Z) in the bottom part corre- 
spond respectively to the experimental phase-matching 
plots in the main and inset panels of Fig. |S| The transi- 
tional wavelengths enumerated above can be easily iden- 
tified there. The SHG ranges where deff 7^ are 1783- 
8224 nm for type-I (ooe) phase-matching in the X-Z 
plane, 2353-6111 nm for type-II (eoe) phase-matching 
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FIG. 10: SHG phase-matching in the principal planes of 
LIS. Thick lines in the lower part show fundamental wave- 
lengths for which doff 7^ and thin lines indicate cases where 
dcff vanishes. The inverse group velocity mismatch GVM 
(A31 — v^^ — v^^ and A32 = v^^ — where Vi,V2,V3 de- 
note the group velocities duoi/dki at Ai, A2, and A3) is shown 
in the upper part for the cases where dcff 7^ 0. The solid 
(dashed) lines correspond to the branch with longer (shorter) 
wavelengths. 

in the X-Y plane, and 2353-2675 nm and 5493-6111 nm 
for type-II (oeo) phase-matching in the Y-Z plane. We 
observe an interesting feature in this crystal: the SHG 
limits with doff ^ are larger for propagation outside 
the principal planes where type-I phase-matching down 
to 1617 nm and up to 8710 nm is possible. Under the 
convention nx < ny < nz the largest birefringence and 
consequently the shortest SHG wavelength is obviously 
achieved for type-I interaction and propagation along the 
Y-axis. However, in LIS when the Y-axis is approached 
in the principal planes X-Y or Y-Z doff for typc-I inter- 
action vanishes and this is true also for the limiting case 
of propagation along the Y-axis. This is the reason why 
propagation outside the principal planes can be used e.g. 
to shorten the SHG lower wavelength limit (Fig. O. A 
similar situation is known to exist and has been experi- 
mentally verified in KTP 69j|. 

It is seen from Fig.^|that in the Y-Z plane the type- 
II interaction is quasi angle-noncritical (but wavelength- 
critical, see inset of Fig. IHl which ensures a large ac- 
ceptance angle and a small walk-off angle. In contrast, 
type-I interaction in the X-Z plane and type-II interac- 
tion in the X-Y plane have regions of quasi-wavelength- 
noncritical phase-matching {d^/dX ^ 0) centered at 
3889 nm and 3803 nm, respectively (Fig. ^1 see also 
Fig. |SJ). Letting $ = ip or 9, the walk-off angles in the 
upper panels of Fig. 1111 are calculated using the simpli- 
fied formula tanpi — [n'i{^)]~^ dnf / valid for uniaxial 
crystals, where the subscript i = 1, 2, 3 is associated with 
wavelength Ai (A;^^ = A^"'^ + A^^) and n^(<I>) is the ex- 
traordinary index of refraction of the walking-off wave as 
given by the uniaxial analogy in the principal planes. We 



FIG. 11: SHG internal angular acceptance (bottom) and 
walk-off angles (top) in the principal planes of LIS. Thick 
solid lines correspond to the branches with longer wavelengths 
from Fig. llOl and thick dashed lines correspond to the branches 
with shorter wavelengths. Only the cases with dcs 7^ are 
included. 

preserved its sign which is in accordance with the corre- 
sponding phase-matching angles i.e. a positive value of 
the walk-off means that the Poynting vector is at an an- 
gle larger than the phase-matching angle and vice-versa. 
Note that the walk-off is very similar for the two branches 
of the solution for phase-matching in all three principal 
planes. In the X-Y plane these branches and the two 
different walk-off parameters (we have here two extraor- 
dinary waves) are almost undistinguishable. The maxi- 
mum walk-off amounts to = -1-1.15° in the X-Z plane 
for SHG at 3889 nm. 

The acceptance angle is evaluated from the phase ve- 
locity mismatch Ak(6^) due to an angular deviation 
(5$ = $ — <I>pM of the wavevectors around the nominal 
phase-matched direction <&pm- A Taylor expansion of Ak 
leads to 

Ak{S<i>) = 7CPM S<i> + 7ncpm(<5*)^ + •■• (35) 

where jcPM—[d{Ak) / 9$]^^^ vanishes for non-critical 
phase-matching (NCPM, i.e™ for $ = 0,90°) and 
7ncpm=^ [^^(^^)/'5'^^]$=$p„- The acceptance an- 
gle, defined as the bandwidth at FWHM of the 
sinc^(AfcL/2) phase-mismatch function, is given by 
A$ — 2.784/|7cpM|i for a critical phase-matching 
(CPM), while for NCPM it is A$ = |2.784/7NCPMi|^/^- 
The lower panels of Fig. ^2 display the acceptance an- 
gles (for a length L = I cm) in the three principal 
planes. The acceptance curves, computed with 7cpm 
only in the expansion H35|l . are interrupted near non- 
criticality at maximum values corresponding to the ones 
calculated for L = 1cm using 7ncpm only. Although 
that type of presentation in Fig. ^] is not very ac- 
curate in these limits we still prefer it since it per- 
mits the results to be presented as scalable with re- 
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FIG. 12; Experimental angular tuning curves for type-II SHG 
of A = 2590 nm in the X-Y plane (at "1>pm = ifiPM = 66.2°, 
blank circles) and in the Y-Z plane ( at $pm = ^pm = 27.9°, 
black circles). The larger acceptance angle in the Y-Z plane 
is due to the much smaller Atiy-z = —0.007 birefringence 
as compared with Anx-Y = —0.035 (Fig. (SJ. The walk-off 
angles are pi ~ ps = 12.2 mrad (X-Y), p2 — 2.8mrad (Y-Z). 
The crystal lengths are L = 6 mm (sample LIS(l) used in 
X-Y) and L — 7mm (sample LIS (2) used in Y-Z plane), see 
Subsection IVIII XI The fundamental beam (diameter 2wo « 
2.6 mm) is the idler output of a ns Nd:YAG-pumped LiNbOs 
OPO. 



spect to the crystal length L (simultaneous considera- 
tion of both derivatives would not allow this). The ex- 
perimental comparison of type-II SHG acceptance angles 
in the Y-Z and X-Y planes is shown with symbols in 
Fig. 113 for A = 2590 nm. The solid lines depict the 
plane-wave sinc^(AfcL/2) phase-mismatch tuning func- 
tions computed with the first-order term in Eq. (|35() , with 
7cpm(6'pm) = (7r/A)[n|p2]epM=7.2 x 10~^fim~^ (type- 
II (oeo) SHG in the Y-Z plane) and 7cpm('Ppm) — 
(VA)[2n^/33-</0i]v=PM= 3.17X 10-Vm-i (type-H (eoe) 
SHG in the X-Y plane). The perfect match of the the- 
oretical curves with the experimental symbols highlights 
both the accuracy of the dispersion data and the single- 
domain feature of the samples used. 

The chosen presentation of the inverse group velocity 
mismatch in the upper part of Fig. llOl is equivalent to the 
spectral acceptance Av but contains the sign as an addi- 
tional information. In the simplest cases of type-I SHG 
or degenerate DFG we have e.g. AvL = 0.886/ IA31I 
and the gain bandwidth in optical parametric amplifi- 
cation assuming a narrow-band pump wave at A3 is in- 
versely proportional to IA21I = \^/v2 ~ The two 
parameters A31 and A32 vanish at 3889 nm in the X- 
Z plane which means large spectral acceptance for the 
SHG of short pulses where the second derivative of the 
wavevector-mismatch comes into play. The situation is 
different for type-II interaction in the X-Y plane: here 
A32 = -A31 at 3803 nm, at 2915 nm A32 vanishes and 
at 4946 nm A31 vanishes but in none of these cases an 



extremum of the spectral acceptance occurs because in 
type-II SHG of short pulses all three waves should be con- 
sidered as broad-band. Thus long crystals can be used to 
increase the conversion efficiency for type-I SHG in the 
X-Z plane near 3889 nm even at femtosecond pulse dura- 
tions. For type-II SHG of femtosecond pulses near 3803 
nm in the X-Y plane the second harmonic will not be 
lengthened either because the temporal walk-off between 
the two polarizations of the fundamental wave will limit 
the interaction length, but consequently the conversion 
efficiency will be low. 



C. Sum-and-difference frequency mixing 

Type-I phase-matching in the X-Z plane for non- 
degenerate three-wave interactions (sum- and difference- 
frequency mixing as well as optical parametric genera- 
tion, amplification and oscillation) is presented in Fig. ?? 
where two branches of the solution can be seen. The 
whole transparency range of LIS can be covered for 
0° < 9 < 40° in the X-Z plane but des increases with the 
phase-matching angle 6. Note that dgff = for NCPM 
(9 = 0°) and that the situation 9 > Vz (Eg.lMjl IS never 
reached. At larger angles 9 we observe a retracing be- 
havior in the left branch: e.g. in the case of OPO, one 
and the same pump wavelength A3 corresponds to two 
pairs (Ai,A2) of signal-idler wavelengths (curve 5). In 
this region the spectral acceptance is very large. Thus 
for 9 = 40° (left branch of curve 5), A3 = 900 nm and 
A2 = 1150 nm all three group velocities are very close: 
at the point where A31 = we have A32 = 23fs/mm. 
This means that this phase-matching configuration is es- 
pecially suitable for frequency conversion of femtosecond 
pulses. In optical parametric amplifiers (OPA's) this ad- 
vantage can be utilized, however, only in combination 
with a control of the spectral bandwidth through the seed 
signal. In the regions near the degeneracy points (SHG 
points) we have on the other hand A21 = but the wave 
at A3 can have in general a different group velocity. Such 
a regime is attractive for broadband parametric ampli- 
fication in the field of a narrow-band pump pulse as in 
the case of chirped pulse optical parametric amplifica- 
tion [t^I. For femtosecond-laser based frequency metrol- 
ogy applications [tH - this phase-matching when imple- 
mented in a synchronously-pumped OPO would provide 
a wide frequency comb grid spanning the 10 — 12 fim 
region useful for high precision spectroscopy of spheri- 
cal reference standard molecules. Increasing the phase- 
matching angle (curve 6 for 9 = 45° in Fig. ??) the two 
branches merge into a closed contour and we approach 
the point where SHG phase-matching only for a single 
wavelength is possible (see Fig. IIOII and all three group 
velocities are again very close but the tunability in that 
case is very limited. 

The curves for type-II phase-matching in the X-Y 
plane (Fig. have a completely different shape. The 
two branches of the solution are represented by curves of 
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FIG. 13: Type-II (eoe and oee) phase-matching for sum- and 
difference-frequency generation in the X-Y plane of LIS and 
several values of the azimuthal angle ip : 35° (lines 1), 40° 
(lines 2), 50° (lines 3), 60° (lines 4) and 90°- NCPM (thick 
lines 5). The curves are terminated at the left and top side 
by the transparency range of the crystal. 



opposite curvature which can cross at two points where 
phase-matching for degenerate DFG or SHG occurs. 
With decreasing phase-matching angle these branches 
separate and a single crossing point is reached (see the 
curves 3 for (p ~ 50°) which corresponds to the single 
SHG solution in Fig. ^| For yet smaller angles no cross- 
ing occurs and the degeneracy point is not reached. For 
all phase-matching curves presented we observe again a 
retracing behavior at longer A3: for each A3 two cou- 
ples (Ai,A2) are phase-matched. At the point where 
these two pairs merge into one (at the maximum A3 per- 
mitting phase-matching) the waves at Ai and A2 have 
equal group velocities and similarly to the case discussed 
for Fig. ?? broadband parametric amplification with a 
narrow-band pump wave can be realized. The deviation 
from the pump group velocity remains, however, essen- 
tial. Thus e.g. at V9 = 90° and A3 = 3310 nm, A21 ~ 
and A31 ~ A32 ~ — 1.65ps/cm. Improved pump group- 
matching at this point occurs when decreasing the phase- 
matching angle ip in accordance with Fig. [TUI 

Comparing Fig.^lto Fig. ?? we note that at relatively 
short A3 (e. g. 1064 nm) full tunability with a fixed crys- 
tal cut is achievable only in the X-Z principal plane. The 
group velocity mismatch depends on the specific wave- 
lengths chosen. At A3 = 1064 nm (a case interesting for 
parametric down-conversion) A2i{X — Z) < A2i{X — Y) 
up to Ai — 10 /im. This means that in this practi- 
cal wavelength range type-II phase-matching in the X- 
Y plane is more advantageous for the development of 
narrow-band parametric generators or oscillators. 

Our analysis of the phase-matching properties of LIS is 
just a first order approximation because it was based on 
the small signal limit where saturation is neglected. We 
considered also only coUinear interactions. Non-coUinear 
interactions attracted recently great interest because in 
addition to the possibility to improve the tunability, to 



spatially separate the beams, and to compensate for the 
spatial walk-off, this can considerably reduce the group 
velocity mismatch |72j . The latter is very promising for 
short pulse interactions. Explicit expressions for non- 
coUinear phase-matching in the principal planes of a bi- 
axial crystal can be found in Ref. |23| and calculations of 
non-coUinear group velocity matching for arbitrary wave 
vector directions in Ref. |68j. 



VIII. SECOND-ORDER NONLINEAR 
COEFFICIENTS 

There has been only one characterization study of the 
nonlinear properties of LIS made by Boyd et al in the 
early seventies 0- The samples used then were of lim- 
ited size and poor quality, and besides their results were 
derived only from non-phase matched SHG at 10.6 /im 
(wedge technique). Further their measurements were 
performed relative to the nonlinearity of GaAs, that was 
several times corrected during the past years. So it is im- 
portant to confirm or to correct these previous results us- 
ing phase-matched second-order interactions on mature 
samples of adequate size. 

The elements of the nonlinear tensor (|27|l have been 
determined from the measurement of the SHG efficiency 
of various laser sources, using either type-I or type-II in- 
teractions in the principal planes of LIS and also out of 
them. Two kinds of measurements have been performed: 
in one kind (Subsection A) , the source is the idler output 
(in the range 2.4 — 2.6 /im) of a 12-ns pulsewidth, 10- 
Hz repetition rate, type-I LiNb03-0P0 19]. Long sam- 
ples were used, and hence the residual absorption and 
the beam walk-off effect, as well as the focusing effects, 
have to be taken into account to derive the nonlinear 
coefficients. The second independent series of SHG ex- 
periments (Subsection B) used femtosecond pulses and 
very thin oriented samples. This allows to neglect the 
effect of absorption and beam walk-off and focusing ef- 
fects, the data analysis in the low-depletion limit being 
derived from the standard plane-wave SHG theory. 



A. Measurements based on Gaussian beam SHG 
theory 

In order to remain in the low pump depletion limit, 
the peak on-axis fundamental intensity used was below 10 
MW/cm^, and the maximum idler energy of the LiNb03- 
OPO was 10 mJ at A = 2590 nm. The spectral linewidth 
of the idler was below 15 GHz, and the transverse beam 
profile ressembled a Gaussian. Three thick samples of dif- 
ferent coloration (rose and yellowish) were used in con- 
junction with weak focusing producing Gaussian beam 
waist of either wq 1mm or wq ~ 1.3 mm. Sample 
LIS(l) was rose-annealed, cut at {ip = 66.5°, 6* ~ 90°) 
for type-II (eoe) SHG in the X-Y plane, with dimensions 
4x4x6 mm^^ (L = 6 mm). Sample LIS (2) was yellowish 
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(as-grown), cut at (</? = 90°, 6* = 28°) for type-II (oeo) 
SHG in the Y-Z plane, with dimensions 4 x 4 x 7mni^ 
(L = 7mm). LIS(l) and LIS(2) were the same sam- 
ples already used in Subsection IVI Al see Fig.|Sland also 
Fig. El Sample LIS(3) had almost identical cut as LIS(l) 
{(fi — 66°, 9 — 90°), but was of yellowish tinge and with 
dimensions 4x4x5 mm'^ (i = 5 mm). We have checked 
that, as the linear dispersion data, the nonlinear coef- 
ficients do not significantly depend on coloration when 
the additional absorption effect of the yellow samples was 
taken into account. For typical pump energy — 1 mJ 
and Wo — 1.3 mm, LIS(2) yielded £'21^ = 2.5 /xJ, while 
LIS(l) yielded E2U, ^llfiJ and LIS (3), = 7^J. The 
difference between the efficiencies of LIS(l) and LIS(3) 
was mainly due to the much stronger second harmonic 
absorption of the latter (a2w = 0.27 cm~^) 0|. For 
LIS(l), the residual fundamental and harmonic absorp- 
tion amounted to a^j 2ui = 0.05 — 0.06 cm"^ while for 
LIS(2) a^,2c. -0.6cm-i [H]. 

The measurements were performed relatively to a 
dual-band anti-reflection-coated KTP sample with L = 
5.93mm, cut at {ip = 0°,e = 56.4°) for SHG of A = 
2.53 /im in the X-Z plane. The measured nonlinear coeffi- 
cient of KTP was ^24 = 2.3±0.2 pm/V J^] and served as 
a reference, without having to correct for dispersion since 
the measurement wavelength is almost identical to that 
used for LIS. We only corrected for the residual absorp- 
tion loss (a^j = 0.068 cm~^ at 2.59 iim) of this KTP sam- 
ple, due to OH^ band W^. According to the Gaussian 
beam SHG theory in the nanosecond pulse regime, the 
effective nonlinearity in the absence of depletion scales 
as Fi 
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where Fsh — -E'2w/-E'5 is the pulse energy conversion effi- 
ciency (in unit of J^^), n is the average index of refrac- 
tion at fundamental and harmonic wavelengths and h is 
the walk-off-focusing function that accounts for diffrac- 
tion, absorption and beam- walk-off effects. The exact 
expression of h for type-II interaction can be found in 
Ref. j76l | . The second harmonic pulse duration was taken 
as T2uj — Ti^/V^. In the nanosecond regime and for 
the given crystal lengths, group- velocity mismatch can 
be safely neglected in SHG. The advantage of measuring 
the LIS coefficient relative to a reference material stems 
from the fact that it is neither necessary to evaluate the 
absolute energies nor to take into account the loss due to 
the numerous filter sets used to isolate the fundamental 
and second harmonic field, since from Ea. (|36|l . the ratio 
of the two defi's scales as 

d2g(LIS) _ rsH(LIS) n2(LIS) L(KTP) /(,(KTP) 



d^giKTP) Fsh (KTP) n2 (KTP) L(LIS) /i(LIS) ' 

(37) 

The efficiencies Fgy were derived from the linear fit of the 
plot of E2UJ versus E'^, for fundamental energies ranging 
from to 5 mJ 0|. The results of the data analysis. 



averaged over the two waist values, are the following: 

doff[LIS(2)] = 1.36 (±15%) deff(KTP), (38) 

deff[LIS(l)] = 3.05 (±15%) doff(KTP), (39) 

deff[LIS(3)] = 3.15 (±15%) dcff(KTP). (40) 

The rather large uncertainties are mainly due to the un- 
certainty in the absorption coefficients through the evalu- 
ation of the h focusing functions. Given that doff (KTP) = 
d24 sine {6 = 58.7° at A = 2590nm), from (gHl and 
(EIJ) one derives rf24(LIS) = 5.73 ± 20% pm/V. The 
two last relations (|39|I - H40|) were averaged to yield, us- 
ing H31|) and the determined value of ^24, (ii5(LIS) = 
7.94 ± 20% pm/V. The ratio of the two nonlinear coef- 
ficients is hence found to be (^24/^15 = 0.72. 



B. Measurements based on femtosecond SHG 

The SHG experiments with femtosecond pulses were 
performed also in the low depletion limit for the funda- 
mental in order to avoid complications from saturation 
effects and spatial effects across the beam cross section. 
The coefficient d^i was estimated by type-I SHG in the 
X-Z plane from Ea. (|33|l . once it was known it was used 
for the determination of d2A from type-II SHG in the X- 
Y plane using Eq. H31|l and assuming Kleinman symmetry 
to hold (i.e. dis = d^i). Finally the diagonal element ^33 
was measured by type-I SHG outside the principal planes 
using Ea. H28l) and the already determined off-diagonal 
components of the d^^^ -tensor, again assuming Klein- 
man symmetry. Note that a similar strategy was used 
previously to determine ^33 of KTP by phase-matched 
SHG 77] and as outlined there the error in the determi- 
nation of this diagonal element is normally larger since 
it accumulates the uncertainties of the other coefficients. 
Alternatively the plane Y-Z could be used for determina- 
tion of (^24 as done in the previous subsection, however, 
the spectral tunability is very narrow in this principal 
plane (see Fig. I10|l which is not adequate for measure- 
ments with broadband femtosecond pulses. 

All three LIS samples used here were annealed and rose 
in colour, had an aperture of 4 x 5 mm^ and were 0.2 mm 
thick. They were cut at {ip ^ 0°,e = 34°) [LIS(4)], 
{if = 59°, 6* = 90°) [LIS(5)], and {ip = 30°, 6* = 44°) 
[LIS (6)], respectively. The small thickness allows to ne- 
glect absorption or scattering losses, the focusing effects 
and also spatial beam walk-off effects. This measurement 
was also relative and we used as a reference sample an 
AGS crystal with the same dimensions cut at (p = 45° 
and 9 = 45° for type-I SHG. In order to minimize the 
error originating from the still different group-velocity 
mismatch the pulsewidth at the fundamental was chosen 
relatively large (160 fs in all cases) since LIS samples with 
thickness less that 0.2mm could not be prepared with suf- 
ficient quality. This pulsewidth corresponds to spectral 
bandwidths of 80-90nm. For dsQ of AGS we used the 
value of 13.9 pm/V jT^ which was measured at a sim- 
ilar wavelength (A = 2.53 /xm) and actually relative to 
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^24 of KTP using the same value for the latter as in the 
previous subsection of the present work. In the plane 
wave approximation and having in mind the equal fun- 
damental energy and crystal thickness, dcs of LIS was 
determined simply from 

dgff(LIS) ^ ^2c.(LIS) [l-fl(AGS)]3 n3(LIS) 
dlg{AGS) £;2a.(AGS) [l-i?(LIS)]3 7i3(AGS) ^ ^ 

where is the product of the the three refractive indices 
involved, and R takes into account the Fresnel reflections 
at the entrance surface (twice for the fundamental) and 
the exit surface (second harmonic). 

The femtosecond source at 1-kHz repetition rate was 
a KTP-based OPA which was seeded by the frequency- 
doubled idler of a BBO-based OPA, both OPAs be- 
ing pumped by the same 800 nm, 40 fs pump source 
(Tirsapphire regenerative amplifler). The 3-mm KTP 
used in the OPA ensured sufRciently long wavelengths 
that cannot be achieved by BBO. Seeding by the fre- 
quency doubled idler of the BBO-OPA was preferred 
against seeding by the signal wavelength for two reasons: 
the signal wavelengths available from the BBO-OPA were 
not short enough to produce idler wavelengths above 
2.5 /im with the KTP-OPA, and frequency doubling in 
a 2-mm thick BBO crystal of the idler used for seeding 
resulted in temporal broadening and spectral narrowing 
which was a prerequisite to improve the accuracy. SHG 
was performed at 2330 nm with LIS (4), 2850 nm with 
LIS(5) and 2300 nm for LIS(6). In all cases more than 
5 /iJ at the fundamental were available so that it was pos- 
sible to reliably measure the SHG with a large area py- 
roelectric detector even at low conversion efficiency. We 
established that up to a conversion efficiency of 20% (en- 
ergy) into the second harmonic the result did not change. 
In the case of LIS(6) the polarization angle 6 was adjusted 
by a He-Ne laser monitoring the polarization rotation by 
the leakage through an analyser till an eigenmode was 
reached. No additional improvement of the SHG effi- 
ciency could be observed later by further alignment of 
this angle which was in good agreement with the calcu- 
lated value from Ea. (|30|l . Note that the dispersion of 
the optic axis angle Vz (see Ea. (|24|l and Fig. (T)) is rela- 
tively small and the deviation between visible (633 nm) 
and 2300 nm is less than 2°, so that the error originating 
from that effect is estimated to be less than 5%. The fol- 
lowing results are values of the d^^^ tensor components 
scaled to A = 2300 nm by using Miller's rule, i.e. assum- 
ing the same wavelength dispersion for the second-order 
nonlinearity as for the linear susceptibility, 

d3i(LIS) = 7.25 (±5%) pm/V, (42) 
d24(LIS) = 5.66 (±10%) pm/V, (43) 
d33(LIS) = -16 (±25%) pm/V. (44) 

The uncertainties are related to the group velocity mis- 
match and to the relatively large angular acceptance us- 
ing thin samples. The above results are in very good 



agreement with the values obtained in the previous sub- 
section by a completely different strategy and we also 
confirmed the different sign of the diagonal element ^33 
as established by Boyd et al Note that if the results 
from Ref. JJ which were obtained by non-phase-matched 
SHG are renormalized using a more recent value of 
di4(GaAs)=83pm/V we arrive at ^31 = 6.14 pm/V, 
'^32 — 5.31 pm/V and ^33 = 9.8 pm/V. Having in mind 
the different wavelength (10.6 /xm) used in Ref. we 
obtain a good agreement within ±20% with the present 
work for all three nonlinear coefficients. 



C. Effective nonlinearity of LIS 

Once the nonlinear coefficients are determined the ef- 
fective nonlinearity dcs can be easily calculated in the 
principal planes using Eqs. I|31|) - (|33|l . The main conclu- 
sion drawn in Ref. [9j, namely that the maximum doff is 
achieved in the X-Y plane where it is wavelength (angle) 
independent because of the similar magnitude of ^24 and 
di5, remains valid with our updated values for the d^ 's. 
However, as already mentioned in subsection lVII Bl inter- 
esting effects can be expected for propagation outside the 
principal planes. Fig. 1141 shows the calculated dcsi'fijS) 
for SHG using Eq. (P5|l - lf^ with the nonlincarities ob- 
tained in the previous subsection and assuming for sim- 
plicity a constant angle Vz = 66° between the optic axes 
and the Z-principal axis (see Fig. U)) which for the wave- 
lengths selected leads to only negligible errors. The fun- 
damental wavelengths indicated in the figure correspond 
exactly to the phase-matching stereographic projections 
presented in Fig.|51where the values for the polar angle 9 
for the limiting cases (p — 0, 90° can be seen. As can be 
seen from the figure for type-I SHG in the Hobden classes 
9, 10, and 11 (see subsection lVII Aj and Ref. 621) the im- 
provement in dcs when propagating outside the principal 
planes is not significant. On the contrary for Hobden 
class 13, it is definitely advantageous to use propagation 
outside the principal planes for maximum d^s in type-I 
SHG. In type-II SHG the curves in Fig.lHhave a differ- 
ent character and maximum dcs is achieved always in the 
X-Y plane which is advantageous as compared to type-II 
SHG in the Y-Z plane. 

For various applications it is important to consider the 
down-conversion into the mid-IR of high power sources, 
e.g. lasers emitting at 1064nm (Fig. I15|l . For the se- 
lected wavelengths, the absolute negative extremum of 
the type-II curves which occurs in X-Y plane (right end, 
with dl'g{90°,Lp) ~ —6.5 pm/V) is larger than the ex- 
tremum of the type-I curves in the X-Z plane (i.e. on the 
left side, with +4 < dX(^> 0°) < 5 Pm/V). Hence type-II 
interaction in the X-Y plane has larger des than type-I 
interaction in the X-Z plane. However, propagation out- 
side the principal planes allows to reach a slightly larger 
maximum —7 < d^^J {9, ip) < —7.5 pm/V for type-I inter- 
action which can be advantageous since other parameters 
like e.g. the spectral bandwidth can be better in type-I 
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FIG. 14: Effective nonlinearity for type-I {ssf, solid lines) 
and type-II {fsf, dashed lines) SHG outside the principal 
planes of LIS as a function of the azimuthal angle ip. The 
labels indicate the fundamental wavelengths corresponding to 
the (ip, 6) phase-matching loci of Fig. |3] 
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FIG. 15: Effective nonlinearity for type-I {ssf, solid lines) 
and type-II (/s/, dashed lines) down- conversion of A3 = 
1064 nm radiation (optical parametric generation, amplifica- 
tion or oscillation) to three selected idler wavelengths Ai = 
4, 7, 10 for propagation outside the principal planes of LIS 
as a function of the azimuthal angle tp. 



interaction. We believe the potential for future applica- 
tions of LIS in the mid-IR is related also to the utilization 
of propagation schemes outside the principal planes. 



IX. PARAMETRIC DOWN-CONVERSION 

A. Femtosecond optical parametric amplification 

As outlined in section LIS possesses the largest 
bandgap among all known mid-IR nonlinear crystals 
which results in an extremely low TPA near 800 nm. This 
unique feature can be utilized for the generation of fem- 
tosecond pulses in the mid-IR by down conversion of the 
most widely used Ti: sapphire femtosecond laser systems 




FIG. 16: Idler spectra (the labels indicate FWHM) that 
demonstrate the achieved tunability with the LIS-based OPA. 
The black circles show the experimentally achieved output 
idler energy, the blank circles show the achieved (internal) 
signal gain, and the diamonds show the applied seed (signal) 
energy. 



operating near 800 nm. In Ref. |21| initial results were 
presented in the 4.8 — 9 /im region but the gain and the 
output energy (2 nJ) were extremely low with a 1.5 mm 
thick LIS sample used as a seeded OPA. Here we report 
an improvement by a factor of about 50 for the output 
idler energy and an extension of the tunability range. 

We applied for this experiment a cube of 5 x 5 x 5 mm"^ 
of rose annealed LIS cut at (/? = 41° and = 90° for 
type-II interaction in the X-Y plane. The single-stage 
travelling-wave-type OPA was pumped by longer pulses 
in order to increase the interaction length. We used 300 fs 
pump pulses at 820 nm from a Ti:sapphire regenerative 
amplifier operating at IkHz. The pulse lengthening was 
achieved by adjustment of the pulse compressor and the 
pump pulses were not bandwidth-limited (appr. three 
times the Fourier limit). About 30 /xJ were used to gen- 
erate continuum in a 2 mm thick sapphire plate and 5 
different interference filters between 875 and 975nm were 
used to select portions (of about lOnm) of the continuum 
for seeding the OPA at the signal wavelength (Fig. I16|l . 
The seed pulse duration was measured by sum-frequency 
cross-correlation with a small fraction of the pump pulses 
in a 0.5 mm thick type-I BBO crystal and deconvolution 
yielded a pulse FWHM of the order of 400 fs (about 3 
times above the Fourier limit). The OPA was pumped 
by 150^ J and the pump beam parameters corresponded 
to a peak on- axis pump intensity of 60 GW/cm^. 

The mid-IR spectra depicted in Fig. 1161 were recorded 
with a 0.5 m monochromator (150gr/mm grating) and 
a liquid-N2-cooled HgCdTe detector. Tunability from 5 
up to almost 12 /im (i.e. the mid-IR cutoff wavelength 
of Fig. [2Jl could be demonstrated. From 8 to 12 ^m an 
almost constant energy level exceeding 80 nJ could be 
achieved, despite the stronger absorption. This can be 
attributed to the lower group velocity mismatch at longer 
idler wavelengths and also to the increasing seed 
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FIG. 17: Cross-correlation traces at two idler wavelengths 
obtained by sum-frequency mixing with a 820 nm, 300 fs ref- 
erence pulse. The deconvolved Gaussian pulse durations 
(FWHM) are in both cases 575 fs. 
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FIG. 18: Wavelength and angle tuning characteristics of type- 
II (eoe) DFG of two cw Ti:Saphire lasers in LIS. The solid 
lines are calculated from the Sellmeier equations of Subsec- 
tion EIZI 



B. Continuous-wave mid-IR DFG 



level which effects compensate for the reduced transmis- 
sion. The pulse duration of the idler pulses was estimated 
from sum-frequency cross-correlation with a 820 nm ref- 
erence pulse in a 0.3 mm thick type-I AGS crystal in 
which the group mismatch between the two input pulses 
ranged from 715fs/mm at Ai = 9.5 iim to 800fs/mm at 
Ai = 6.5 /xm. Fig. El shows the corresponding cross- 
correlation traces. The deconvolved pulse duration at 
7.65 and 9.51 /im (585 fs in both cases) corresponds to 
bandwidth-limited Gaussian shaped idler pulses (time 
bandwidth product equal to 0.44). 

Thus LIS is the only nonlinear crystal that permits 
the direct frequency conversion of short pulses from the 
800 nm region to the mid-IR above ~ 5 fim in a sin- 
gle step extending the region covered by crystals like 
LiNbOs, KNbOa or KTP and its isomorphs. The re- 
sults in this section were obtained with amplified pump 
pulses at low (IkHz) repetition rates. Similar single 
stage down conversion of high repetition rate (100 MHz) 
Tiisapphire femtosecond oscillators should be possible in 
synchronously pumped LIS-based OPOs. 

A further potential application of LIS in femtosecond 
technology which relies on the extremely low TPA near 
800nm [21| is related to the generation of single (or sub-) 
cycle femtosecond pulses (both at low and high repetition 
rates) in the mid-IR by DFG of spectral components be- 
longing to the same broadband spectrum of an ultrashort 
800nm pulse f^. If 15 fs long bandwidth-limited pulses 
centered at 800nm are used for this purpose the generated 
mid-IR wavelength is about 10 /xm at which wavelength a 
30fs long pulse would possess only one optical cycle. The 
basic limitation in the conversion efficiency when u sing 
GaSe as previously demonstrated for this purpose |79| 
was its strong TPA at SOOnm. 



We report in this subsection the first DFG experiments 
(5.5 < Ai < 11.3 /zm) ever performed with LIS, whose 
partial results were already used to derive the accurate 
Sellmeier equations of Subsection lVI Al Continuous- wave 
(cw) DFG is most suited for high-resolution spectroscopy 
in the mid-IR range: when both DFG laser sources have 
linewidths < 1 MHz the mid-IR radiation at Ai allows to 
perform sub-Doppler spectroscopy of heavy atmospheric 
molecules with only ~ 1 fj.W of power. Until now, cw 
DFG of two near-IR lasers for the generation of deep 
mid-IR fight either used AGS |80l (sj, |82| or GaSe [Hi . 
We demonstrate that LIS is an alternative promising can- 
didate, since we report cw parametric generation up to 
11.3 /im with an uncoated yellowish LIS sample of dimen- 
sion 5x5x10 mm"^ (L = 10mm), cut in the X-Y plane for 
type-H (eoe) DFG at {e = 90°,ip = 42°). 

Two tunable cw single-frequency Ti:sapphire lasers 
were used as difference-frequency mixing sources. The 
laser beams, orthogonally polarized, were collinearly fo- 
cused onto the crystal with a 35-cm-focal length lens, 
yielding a nearly optimal wq — 55 /im waist for both 
pump wavelengths A3 < A2. The pump wavelengths 
were blocked after the sample by a 1-mm thick uncoated 
germanium (Ge) filter and the mid-IR radiation was de- 
tected by a calibrated liquid-N2-cooled HgCdTe photo- 
conductive detector with a 1 x 1 mm^ active area. The 
detected voltage was amplified by a low-noise preampli- 
fier and fed to a lock- in amplifier. 

In Table IIVI we have already seen that at fixed 
(p = 42°, wavelength tuning alone can cover only the 
6.6 — 7/im range, limited by the tunability range of the 
Ti:sapphire lasers (from curves 2-3 of Fig. ^1 a tunability 
range extending into the blue would be required to ex- 
tend the DFG to ~ 10 nm). Taking profit of angle tuning, 
mid-IR radiation was generated with broad wavelength 
tunability in the 5.5 — 11.3 /zm spectral region with com- 
bined pump laser wavelength (A3) and crystal angle (ip) 
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FIG. 19: Mid-IR DFG spectral tuning curve in pump tuning 
mode around A3 ~ 765 nm. The other wavelength was fixed 
at A2 = 865.3 nm. In this frequency scan, we could observe 
atmospheric water vapor absorption lines in laboratory ambi- 
ent air over 42.5-cm open path between the LIS crystal and 
the HgCdTe detector. 



tuning [sj: in Fig. ll8lone Ti: sapphire laser was fixed at 
A2 ~ 836 nm, and the other (A3) was tuned from 724 to 
778 nm. The sohd lines are calculated phase-matching 
curves from Eas. H21|) - (|23|l . highlighting the accuracy of 
the Sellmeier equations that were fitted using the fixed- 
angle DFG data of Table Hvl 

Fig. 1191 displays the mid-IR spectral tuning bandwidth 
around 1510 cm~^ (Ai ^ 6.6/im), when A2 is set at 
865.3 nm and A3 is tuned around 765 nm, corresponding 
to a quasi-normal incidence phase-matching. A APi ^ 
lOcm^i {v = l/X = vjc) FWHM spectral bandwidth is 
obtained, within which sharp atmospheric water vapor 
absorption lines can be observed. The FWHM plane- 
wave spectral acceptance for this pump tuning mode ex- 
presses as 



l\vx = 1.722- 
L 



dAk 



(45) 



where dAk/dvi can be expanded, taking the constraint 
1^3 ^ 1^2 ~ 1^1 into account, as 



dAk 27r 



ns,{if>) - A3 



ni{(p) + A 



dXi 



27rA 



31- 



(46) 



In Ea. (|46|l . the angle-dependent indices are the extraor- 
dinary indices in the X-Y plane. From (|45|I - H46|I . the 
calculated spectral acceptance is Ai/i = 214 GHz (A^i ~ 
7.2 cm~^) corresponding to a pump tuning over AA3 — 
0.42 nm. The bandwidth broadening observed in Fig. 1191 
is attributed to the focused beam effect . 

Let us note that despite the strong mid-IR absorption 
at 11.3 /im (Fig.EJl, we could still detect 3nW mid-IR 
power for input powers P3 = 200 mW and P2 — 265 mW. 



This power is generated by the last exit layers of material. 
For the same incident powers, the detected power at Ai = 
7/im was Pi = 33.5 nW. The net conversion efficiency 
at this wavelength, corrected for all Fresnel optical loss, 
amounted to F = P1/P3P2 ^ 1.6fiW/W^. Using the 
general focused-beam theory for DFG 84] that evaluates 
the focusing function h for type-II DFG for the given 
focusing parameters, walk-off value (|pi,3| — 0.97°) and 
absorption loss (03 ~ 0.15 cm~^, a2 ^ 0.1 cm~^, ai ^ 
0.05 cm"""^), we derive an effective nonlinear coefficient 
dcs(J fim) ~ 5.3pm/V. This value is slightly lower than 
the calculated value using Ea. (|^ and the determined 
values of ^24 and di5 (Eqs.||32|)-(|llll), dcs ~ 6.55 pm/V. 
The difference may arise from the residual atmospheric 
water absorption lines in this region. 



X. 



LASER-INDUCED DAMAGE IN LilnSo 



Illumination with various laser sources resulted in the 
following photo-induced phenomena observed in LIS, 
whose magnitude depends on wavelength and pulse du- 
ration: (i) micro-cracking (at A = 5 ^m using femtosec- 
ond pulses from a free-electron laser [l^); (ii) dendrite 
structure formation (at A = 1.079 /zm using nanosecond 
pulses); (iii) photo- induced absorption (PI A) as grey- 
tracks (at A = 0.8 /im using femtosecond pulses). Op- 
tical damage of (i) and (ii) kinds is irreversible, while 
PIA damage is to a certain extent reversible. It disap- 
pears at lower levels spontaneously after several hours 
of storage at room temperature or several minutes of il- 
lumation with A > 500 nm. The PIA was detected first 
in femtosecond CPA experiments with a 1.5mm thick 
LIS sample pumped by a IkHz Ti:sapphire regenerative 
amplifier delivering pulses at 800nm with a duration of 
200fs and energy of 200 J . No surface damage was 
observed in these experiments up to a peak on-axis inten- 
sity of 140 GW/cm^. In principle as with other materials 
optical damage in the femtosecond regime is expected 
to occur indirectly as a result of filamentation and self- 
focussing through higher-order nonlinear effects. These 
studies permitted to estimate only an upper limit on the 
nonlinear absorption losses at 800 nm which were inde- 
pendent of polarization and presumably caused by inclu- 
sions. They could not be fitted by a TPA process and 
the estimated upper limit of TPA is /?tpa ~ 0.04cm/GW 
at 800 nm 2]^. This is about 100 times lower than the 
Ptpa = 3.5 — 5.6cm/GW (depending on the polariza- 
tion) fitted for a 0.3mm thick AGS crystal under identi- 
cal conditions, which is a manifestation of the advantages 
of LIS related to its larger bandgap. An updated value 
for the upper TPA limit that we derived for a different 
sample, namely the LIS crystal used in Subsection IIX Al 
is only slightly higher: /3tpa — 0.05cm/GW. Note that 
recent TPA measurements on KNb03 (which belongs to 
the same class mm2) at 846 nm using a 5-mm thick sam- 
ple phase-matched for cw blue light generation yielded a 
TPA magnitude (/3tpa = 3.2cm/GW) of the same order 
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as that for AGS [83 ■ In the experiment of Ref. zl] the 
gray tracks formation (leading to PIA) was attributed 
to non-phase-matched SHG and residual blue light ab- 
sorption, i.e. this damage mechanism is only indirectly 
related to the high intensity of the femtosecond pulses 
that resulted in higher probability for unwanted SHG. 
Subsequently differential absorption (before and after il- 
lumination) was measured that had a magnitude reaching 
several cm~^ and extended from 400 nm up to the near- 
IR (1 — 2 /j,m) 01 . Well pronounced maxima of this ab- 
sorption were observed at 430 nm and 650 nm. This pho- 
tochromic induced-absorption is similar to the so-called 
Blue-Induced Infrared Absorption (BLIIRA) observed in 
SHG experiments with KNbOa, for which an increased 
near-IR absorption coefficient (0.011 to 0.056 cm~^ de- 
pending on the sample and corresponding to 2.5 to 9 
times the usual linear absorption) at the fundamental 
wavelength of 860 nm was detected with 100 mW of blue 
SH power (30kW/cm^) generated with a Ti:sapphire 
laser [s^j- The BLIIRA damage was also reported to be 
reversible, after a few hours of relaxation at room tem- 
perature, and could be minimized by heating the crystal 
or pumping at longer wavelengths [88|. 

In order to understand the origin of the PIA in LIS, 
we proceeded with a direct illumination of LIS samples 
at 300 < A < 450 nm using a 1 kW Xe lamp filtered 
through a monochromator for ~ 10 minutes: The same 
coloration effect was indeed obtained, with a maximum 
effect at 360 nm excitation and vanishing at wavelengths 
shorter than 300 nm. Thus the PIA observed in the OPA 
experiments is really caused by linear absorption of non- 
phase matched blue SHG by point defects present in the 
lattice. In line with the photo-induced signal in absorp- 
tion spectra (see Fig. 3 in Ref. 86]), a complicated system 
of lines appears in the Electron Spin Resonance (ESR) 
spectrum. Analysis of the angular dependence allowed 
us to reveal the superposition of two spectra related to 
different paramagnetic defects. One of them has a spin 

= 1/2 and an anisotropic g-factor with the following 
principal values of the g-tensor: gi — 2.160, 52 = 2.023, 
gs = 2.023. This spectrum is associated with a center 
of axial symmetry with gi, directed along the In-S bond 
and is related to sulfur S~ ion in a regular site. The sec- 
ond spectrum comprising a broad line at if ~ 3300 Gauss 
{g = 2.0024) is associated with Li°. 

The interpretation of these results is the following. At 
short-wavelength excitation (A < 450 nm) the following 
recharge process occurs in LIS: 



Li^ 



Li^ 



(47) 



Such excitation produces free charge carriers which are 
captured in some deep traps and can be released at 
A > 500 nm illumination. We established that the gray 
track formation is maximized for excitation near 360 nm 
but dropped by one order of magnitude at 420 nm, which 
means that in OPO/OPA experiments it can be reduced 
by using slightly longer pump wavelength A3. Note that 
the same process also takes place at x-ray irradiation, 
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FIG. 20: Measured pulse energy with crystal inserted 
(squares) and without it (diamonds). The drop in transmit- 
ted energy (last square) indicates damage formation. Inset: 
Damage threshold, expressed as fluence, versus beam radius 
at the crystal surface. 



which is used for crystal orientation when fabricating op- 
tical elements: some darkening occurs sometimes in the 
bulk sample during this operation. The PIA magnitude 
was found to be dependent on the LIS composition j89j : 
The effect is minimum in the case of nearly stoichiometric 
composition. The PIA was also found to be much more 
pronounced in LIS than in LISe. 

Systematic optical damage threshold measurements 
were performed with Nd:YAG lasers at 1064 nm - i.e. 
well within the clear transmission range and far from the 
TPA region. Initial measurements with a 6 W cw focused 
Nd:YAG laser beam failed to produce measurable dam- 
age, and subsequent measurements were performed with 
a pulsed laser. The pulse energy was varied between 6 
mJ and 25 mJ at a repetition rate of 10 Hz and a pulse 
width of about 10 ns. The coUimated Gaussian beam 
with a beam radius of w,;„ = 1.47 mm was focused with a 
200 mm focal length lens to produce a beam waist with 
a radius of wo= 46 /im. The crystal used was a rose sam- 
ple with dimensions 4x5x5 mm'' {9 = 90°,(p ^ 24° ) . The 
corresponding fluences ranged from 0.19 J/cm^ to 752 
J/cm^, and by using different pulse energies, the same 
fluence could be sustained over a range of spot radii. 
The Gaussian character of the beam was verified, and 
the beam radius measured by the knife-edge technique. 
The onset of radiation damage was determined from the 
observed sudden drop in the transmitted signal as shown 
in Fig. 1201 At each experimental point, the crystal was 
inspected visually in a low-power microscope, and the 
drop in transmission was found to correlate with the on- 
set of visible damage at the surface. By performing ex- 
periments with different pulse energies it was found that 
the damage threshold, expressed in J/cm^, tends to de- 
crease with increasing spot size as seen in the inset figure 
of Fig. [5D| Although the amount of data is too limited 
to warrant a detailed analysis, it suggests a levelling off 
at about 1 J/cm^ at large spot sizes (i.e. a peak inten- 
sity threshold of 100 MW/cm^), and this value agrees 
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well with the value 1.1 J/cm^ reported by Knippels et al. 
at a wavelength of 5 /im |23|. The observed dependence 
in the inset figure is understandable in view of the fact 
that as the spot size increases, it becomes increasingly 
difficult for heat generated at the centre to diffuse to the 
unexposed regions of the crystal. In Fig. [201 the initial 
increase of transmitted energy reflects the fact that the 
center of the beam is too close to the edge to allow the 
entire beam to be intercepted by the crystal when the 
beam radius is largest. In some cases damage was found 
to extend into the interior of the crystal. 

These damage properties should be discussed in com- 
parison with available damage data on AGS or AGSe. 
For AGS under cw pumping the onset of thermal lensing 
has been previously evaluated at ^ llkW/cm^ from a 
cw AGS doubly resonant OPO containing ^ 3 W of cir- 
culating power at 1.26 and 2.52 /im focused to a waist 
of wq = 90 /^m 1^ 0] . The cw damage threshold value 
is actually unknown, but is estimated to occur in the 
range - 20 - 40kW/cm2 despite a report that 3.5 
W of cw Nd:YAG radiation focused to a waist of 20 /im 
(280kW/cm2) did not damage an AGS sample For 
AGSe, the measured damage threshold at 9 um (cw CO2 
laser) was ascertained to be 5 — 22 kW/cm^ [221 > corrobo- 
rating that chalcopyrites display extremely low cw dam- 
age threshold. The fact that for 6W at 1064 nm, focused 
down to Wq — 40 /im no visible damage was observed in 
the LIS sample sets the cw damage threshold to a lower 
limit of > 120kW/cm^. The 5 times larger thermal con- 
ductivity of LIS as compared to AGS certainly favors the 
use of LIS in intra-cavity cw applications. Under similar 
nanosecond pulsed regime, however, the ~ 1 J/cm^ dam- 
age threshold measured with LIS is c omp arable to the 
threshold found for AGS or AGSe [H IJ . Independent 
measurements at longer wavelengths (9.55 /xm) with 30 ns 
long pulses yielded, under identical experimental condi- 
tions, damage thresholds of 6.5J/cm^ for LIS jH^I and 
5.5J/cm^ and 4.2 J/cm^ for AGS and AGSe, respec- 
tively m. 

From these preliminary damage threshold investiga- 
tions, it is premature to assess whether the damage mech- 
anisms are due to native point defects or are intrinsic to 
the compound. Considering the current state-of-the-art 
of residual absorption on this compound, one may ex- 
pect that with an improvement in residual transparency 
and a better control of the growth parameters toward the 
production of stoichiometric material, the above damage 
threshold values can be improved. 



XI. CONCLUSIONS 

We have characterized the main physical properties 
of the newly emerging lithium thioindate chalcogenide 
semiconductor that are relevant for nonlinear applica- 
tions from the visible up to the deep mid-IR. Table IVII 
summarizes these properties, as well as those determined 
in previous works. LiInS2 or LIS can now be consid- 



ered as a mature compound and its future widespread 
use in nonlinear devices should refine all the data col- 
lected from this investigation campaign. Among its most 
relevant advantages over the existing mid-IR crystals, one 
may quote its excellent thermal stability and high dam- 
age threshold suitable for high power laser applications 
and its extremely extended phase-matching capabilities 
(over its whole transparency range) using either type-I 
or type-II interactions in or out-of principal planes. Such 
an extended phase-matching capability is still compat- 
ible with low birefringence walk-off angles. LIS is the 
only existing material for direct down-conversion from 
the near-IR to the deep mid-IR up to 12 /xm, owing to 
its highest bandgap energy. Its effective nonlinearity of 
~ 7pm/V (slightly lower than that of the chalcopyrite 
crystals) is substantially higher than that of the oxides 
of the KTP family and is almost constant for type-II 
interaction in the azimuthal principal plane. We have 
shown that LIS should be very promising for the fre- 
quency down-conversion of femtosecond laser sources. At 
equivalent loss level, LIS would be less subject to thermal 
lensing than AGS 0, [3 permitting its use in cw OPOs 
for deep mid-IR generation. 

The next challenge and effort will be directed to the 
detailed studies on its eventual ferro-electric properties. 
In contrast with the chalcopyrites, LIS and LISe are py- 
roelectric as the orthorhombic symmetry oxides (KTP 
isotypes) belonging to the same mm2 point group. The 
spontaneous polarization is due to the fact that in each 
bipyramid formed by two adjacent tetrapores (Fig. 
only one is occupied by a Li+ cation. The other one is 
empty, which corresponds to a negative effective charge. 
An intriguing issue, which has not still received a definite 
answer, is whether such Li-chalcogenides can be related 
to ferroelectric materials, i.e. if it is possible to invert 
their spontaneous polarization along the polar c-axis. In 
their short note |53| , Negran et al concluded that LIS and 
LISe do not display a para-to-ferroelectric phase transi- 
tion from the absence of a sudden anomalous rise in the 
record of the pyro-electric current versus temperature. 
This is not surprising because of the screening effect of 
the high ionic conductivity of Li. They measured a pyro- 
electric coefficient p'^(300X) = 2.6 x 10-i°C/(cm2.K). A 
more recent and reliable measurement of the pyro-electric 
coefficient performed with our rose samples yielded a 2.3 
times higher value f Table jVI)) . which is 5 times smaller 
than the value for KTP _96J. In contrast with Negran et 
al 's report, a strong increase of the pyro-current slightly 
above room-temperature, preceded by a sharp current 
spike, was noticed but this increase was attributed to 
space c harg e relaxation of the thermally-activated Li+ 
cations |55| . 

The issue about whether Li-chalcogenides are (or are 
not) related to ferroelectrics is of upmost importance, 
considering their deep mid-IR transparency which would 
allow to fabricate quasi-phase matched periodically-poled 
structures for wavelengths longer than ~ 4 — 5 /^m, which 
is the current limitation of oxide-based ferroelectrics such 
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TABLE VI; Summary of the main physical properties of LiInS2. 



CRYSTALLOGRAPHIC DATA 
Structure 

Symmetry, point group 
Space group 
Lattice parameters (A) 
Principal axes assignment 
Density (g/cm'^) 
Microhardness (Pa) 



OPTICAL PROPERTIES 
Optical transmission d-tm) 
Bandgap energy (T = 300 K) (eV) 
Indices of refraction at T = 20° : 

A = 0.532/im 

A = 1.064/im 

A = 10.6^m 
Optical axis angles (0.35/im< A < 11.5/im) 
Birefringence walk-off at 1064 nm (mrad) 
SHG fundamental wavelength range (nm) 
Thermo-optic coefficients (10~^/°C) 

at 1064 nm, T = 20°C 
Absorption coefficient at 1.064 /im (cm^^) 
Laser damage threshold: 

• surface, at 10 ns, / = 10 Hz, A — 1.064 /im 

• bulk, at 500 fs, / = 25 MHz, A = 1.064 /xm 

• surface, at 200 fs, / = 1 kHz, A = 0.8 ^m 
Two-photon absorption at 0.8 (cm/GW) 
Nonlinear optical coefficients (pm/V) 

• by cw SHG at 2.5 

• by femtosecond SHG at 2.3 /xm 



Wurtzite-type (/J-NaFeOa structure) 
Orthorhombic, mm2 (negative biaxial) 
Pna2i 

a = 6.874; b = 8.033; c = 6.462 (Ref. [18]) 

{X,Y,Z) ^ {b,a,c) 

3.46 

145x10'^ (Refs. [26. 27]) 



0.34 - 11.5 

3.55 (E II c); 3.61 {E \\ b); 3.58 {E \\ a) (Ref. [18]) 

nx = 2.2413; ny = 2.2824; m = 2.2991 

nx = 2.1288; ny = 2.1645; m = 2.1734 

nx = 2.0200; ny = 2.0566; m = 2.0595 

55° <Vz < 80° 

px-Y = 14.7; px-z = 20 

1617 - 8710 

dnx/dT = 3.72; dny/dT = 4.55; dnz/dT = 4.47 
a < 0.04 

1 J/cm^ 

> 6 GW/cm^ (Ref. [20j) 

> 140 GW/cm^ (Ref. [21j) 
/3tpa < 0.05 

di5 = 7.9(±1.6); d24 = 5.7(±1.2) 

dai = 7.2(±0.4); ^24 = 5.7(±0.6); ^33 = -16(±4) 



THERMAL PROPERTIES 
Melting point (°C) 

Thermal expansion at T = 20°C (10"V°C) 
Specific heat at T = 300 K (J/mol/K) 
Thermal conductivity at T = 300K (W/m/°C) 



ELECTRICAL PROPERTIES 

Pyroelectric coeflicient at T = 300 K {pC/m^K) 

Pure electro-optic coefficients (pm/V) 

Piezo-electric coefficients (pm/V) 

Electric conductivity (f2~^cm^^) 



r„,eit ~ 1000 

ax = +1.64; ay = +0.91; qz = +0.68 
C = 92.9 

Kx = 6.2; Ky = 6.0; i^z = 7.6 (Ref. [50]) 



6 (Ref. [55j) 

ri3 = 1.0rr23_= 0.4; r33 = -1.3 
d3i = -5.6; d32 = -2.7; ^33 = +7.5 

az = 4.3 X 10"^° (Ref. [98]); cr = 1 x 10"^^ (Refs. [22,^) 



as LiNbOs or KTiOPOs. The measurement of the Curie 
temperature (a temperature at which the material tends 
to loose its spontaneous polarization), the study of its 
dielectric conductivity versus an applied ac electric field 
in order to reveal a dielectric hysteresis loop (and the 
value of the coercive field) have to be performed to answer 
to the question. 

Finally, we started to investigate other Li-containing 
compounds such as LiInSe2 (LISeL LiGaS2 (LGS), 
LiGaSea (LGSe), LiGaTci (LGTe) 
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